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ABSTRACT 

Mineralogical studies of silicate features emitted by dust grains in protoplanetary disks and Solar 
System bodies can shed light on the progress of planet formation. The significant fraction of crystalline 
material in comets, chondritic meteorites and interplanetary dust particles indicates a modification 
of the almost completely amorphous interstellar medium (ISM) dust from which they formed. The 
production of crystalline silicates, thus, must happen in protoplanetary disks, where dust evolves to 
build planets and planetesimals. Different scenarios have been proposed, but it is still unclear how and 
when this happens. This paper presents dust grain mineralogy (composition, crystallinity and grain 
size distribution) of a complete sample of protoplanetary disks in the young Serpens cluster. These 
results are compared to those in the young Taurus region and to sources that have retained their 
protoplanetary disks in the older Upper Scorpius and 77 Chamaeleontis stellar clusters, using the same 
analysis technique for all samples. This comparison allows an investigation of the grain mineralogy 
evolution with time for a total sample of 139 disks. The mean cluster age and disk fraction are used 
as indicators of the evolutionary stage of the different populations. Our results show that the disks 
in the different regions have similar distributions of mean grain sizes and crystallinity fractions (~10 
- 20%) despite the spread in mean ages. Furthermore, there is no evidence of preferential grain sizes 
for any given disk geometry, nor for the mean cluster crystallinity fraction to increase with mean age 
in the 1-8 Myr range. The main implication is that a modest level of crystallinity is established 
in the disk surface early on (< 1 Myr), reaching an equilibrium that is independent of what may be 
happening in the disk midplane. These results are discussed in the context of planet formation, in 
comparison with mineralogical results from small bodies in our own Solar System. 

Subject headings: stars: pre-main-sequence ~ planetary systems: protoplanetary disks - circumstellar 
matter - infrared: stars - methods: statistical 



1. INTRODUCTION 

Protoplanetary disks originate from dense cloud ma- 
terial consisting of sub-/zm sized, almost completely 
amorphous interstellar medium (ISM ) dust grains 

; (iBeckwith et al.l [20001 : Hi fc Draindl25oit iKemper et al.l 

. 120041: lHenningll2010f ). The dust and gas in these disks 
form the basic matter from which planets may form. 

\ At the same time, mineralogical studies of primitive so- 
lar system bodies suggest that a considerable fraction of 
the s ilicate grains in these objects are of crystalline na- 

; ture ([Wooden et al.ll2007tlPontoppidan fc Brearievll2010L 
and references therein) . It is then naturally implied that 
the crystallinity fraction increases, through thermal and 
chemical modification of these solids during the general 
planet formation process, commonly referred to in the 
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literature as "disk evolution" . 

As time passes, the small dust responsible for the in- 
frared (IR) excess observed around young stars is sub- 
jected to different processes that affect, and will eventu- 
ally determine how this progression will end. Planets and 
planetary systems have been observed around hundreds 
of stars other th an the Sun, showing that this result is 
rather common (|Udry fc Santo£|[2007D . IR observations 
have revealed a great number of debris disks, composed 
of large planetesimal rocks and smaller bodies, around a 
variety of stars spa nning a large range in spectral types 
and a ges (iRieke et al . 2005; Brvden et al. 2006; Su et al] 
[2OOI ICautier et al.ll2007l : (Carpenter et al.1l2009t ). A few 
debris disks are kno wn to harbor planets (e.g., Pic- 
toris and Fomalhaut, ILagrange et al.l [20101 : IKalas et al.l 
|2008[ ). although it is still unclear whether this is often 
true ([Kospal ct al. 2009). The majority of main-sequence 
stars show no signs of planets or debris within the cur- 
rent observational limitations, however, indicating that 
the disks around such stars at the time of their forma- 
tion have dissipated completely, leaving no dust behind 
to tell the story. Which processes are important and 
determinant for the aftermath of disk evolution are still 
under debate, and this topic is the subject of many the- 
oretical and observational studies over the last decade, 
stimulated in large by recent IR and (sub-)millimeter fa- 
cilities. 

Specifically on the subject of the mineralogical compo- 
sition, spectra from the ground and the Infrared Space 
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Observatory gave the first clues of a potential link be- 
tween crystalline material in protoplanetary disks and 
comets. A great similarity was noted between the spectra 
of the disk around t he Her big star HP 100546 and that 
of comet Hale-Bopp (jCrovisier et al.iri997l : IMalfait et alj 
[T9981) . M ore recently, the In fraRed Spectrograph (IRS, 5 
- 38 fim, iHouck et al.ll2004 ) on-board the Spitzer Space 
Telescope allowed an unprecedented combination of high 
sensitivity and the ability to observe large numbers of 
disks, down to the brown dwarf limit. The shape of the 
silicate features probed by the IRS spectra at 10 and 
20 /xm is affected by the composition, size and struc- 
ture of its emitting dust. Amorphous silicates show 
broad smooth mid-IR features, while the opacities of 
crystalline grains show sharp features due to their large- 
scale lattice arrangement, such that even small fractions 
of crystalline gra ins produc e additio nal structure in the 
silicate features (IMin et all [2005.: Bouwman et al.l [20081 : 



lJuhasz et al.l 12009: Olo fsson et al.il2oToi r Because most 
protoplanetary disks are optically thick at optical and IR 
wavelengths, the silicate features observed in the mid- 
IR are generally emitted by dust in the optically thin 
disk surface only. To probe the disk midplane, observa- 
tions at longer wavelengths are necessary. Additionally, 
the emission at 10 and 20 fim has been shown to arise 
from different grain populations, probing different rad ii 
(|Kessler-Silacci et al.l [20061: lOlofsson eran[200l [20T0I) . 
While the 10 /xm feature probes a warmer dust popu- 
lation, at < 1 AU for T Tauri stars, the dust emitting 
at 20 nui is colder, further out and deeper into the disk 
(jKessler-Silacci et al.|[2007[) . 

Two methods have been proposed to explain the for- 
mation of crystal grains: thermal annealing of amor- 
phous grains or vaporization followed by gas-phase con- 
densation. Both methods require high tempe ratures 
(above -1000 K, iFabian et all l2000t [Gail [200I which 
is inconsistent with outer disk temperatures. However, 
crystalline grains have been observed in outer, as well 
as in inner disks (jvan Boekel et al1l2004[ ). Large-scale 
radial mixing has been invoked to explain the pres- 
ence of crystals at low temperatures i n the outer disk 
(|Bockelee-Morvan et al.i i2000: Gail200l[CiesIa.2009 '). A 
third proposed formation mechanism for crystal forma- 
tion is that shock waves could locally heat amorphous 
silicates and crystalli ze them ([Desch fc ConnoUvl 120021 : 
IHarker fc Desch'2002^. 

Prom protoplanetary disks to comets, several authors 
have attempted to infer the dust composition from IRS 
spectra and laboratory data on amorphous and crys- 
talline silicate dust, using a variety of analysis tech- 
nique s. Whether for individual objects (Forres t et al. | 
20041: IMerfn et all ,2007; Pintc ct al. 2008; Bouv et al.l 



20081). for mixe d disk samples (Bouwman ct al. 2001'; 



Apai et aLll2005l: [van Boekel et al. 2005: Bouwman et al 



20081: lOlofsson et al.ll2009i r2010l: iJuMsz et al.ll2010ir or 

systematic studies of the di sk population of a given star- 
forming reg ion ( Sicilia- Agui lar et al.ll2009l : iWatson et al.l 
120091 : [Sargent et al.ll2009 ). it has been shown that a sig- 
nificant mass fraction of the dust in those disks must be 
in crystalline form. However, the many studies dealing 
with the mineralogical composition of dust to date fo- 
cus on a specific region or object, failing to investigate 
the hypothesis that the crystallinity fraction is a mea- 
sure of the evolutionary stage of a region. That is, no 



study in the literature has yet investigated an increase of 
crystallinity fraction with cluster age. 

Mineralogical studies of Solar System bodies show a 
range of crystallinity fractions. Evidence from primitive 
chondrites shows that the abundance of crystalline sili- 
cate material varies from nearl y nothing up to 20 - 30 % 
(e.g. Acfer 094 and ALH77307. [Pontoppidan fc BrearlevI 
I2OIOI and references therein). Oort cloud comets, with 
long periods and large distances from the Sun, have 
inferred crystallinity fractions up to 60 - 80 % (e.g. 
Hale-Bopp, Wooden _et_ al. 1999, 2007j). Jupiter-family, 
or short period comets, have lower fra.ctions , up to —35 
% (e.g. 9P/Tempel 1 IHarker et al.l [2007t 81P/Wild 
2, iZolenskv et al.ll2006l) . This discrepancy in fractions 
points to the existence of a radial dependence in crys- 
tallinity fra ction in the protop lanetary disk around the 
young Sun ([Barker et al.|[2005[ ). It is important to note 
that those values are model dependent, and the use of 
large amorphous grains (10-100 ^ni ) can lead to systen i- 
atically lower crystalline fractions (IHarker et al."2002f). 
This is evident for Hale-Bopp, where IMin et al. (20051) 
find a much lower fraction (—7.5 %) than other authors, 
using a distribution of amorphous grain sizes up to 100 
fim. What is clear is that even within the discrepancies, 
the crystallinity fractions derived for Solar System bod- 
ies are appr eciably higher than t hose derived for the ISM 
dust (< 2%. lKemper et al.]|2004[ ). Recent Spitzer data in- 
dicate further similarities between crystalline silicate fea- 
tures seen in comets or asteroids with those seen i n some 
debri s disks around solar m ass stars (Beichman et al.l 
2006: Lisse et al.ll2007l [20081) . One proposed explanation 
is that the observed spectral features in the disk result 
from the catastrophic break-up of a single large body 
(a 'super comet') which creates the small dust particles 
needed for detection. At the even earlier protoplane- 
tary disk stage, there is limited observational evidence 
for radial gradients in crystallinity from mid-infrared 
interferometry data, with highe r crystallinity fra,ctions 
found closer to the young stars ("van Bo ekel et all 120041 : 
Schcgcrcr ct al. 2008). All of this suggests that the crys- 
tallization occurs early in the disk evolution and is then 
incorporated into larger solid bodies. 

Besides dust composition, the evolution of grain sizes 
is an essential indicator of disk evolution. The initially 
sub-/im size ISM grains must grow astounding 14-15 or- 
ders of magnitude in diameter if they are to form plan- 
ets. If grains were to grow orderly and steadily, theo- 
retical calculations predict disks to have fully dissipated 
their small grains within —10^ y ears (jWeidenschillind 
[T980HDullemond fc Dominia[200g) . The fact that many 
disks a few Myr old are observed to have small grains 
(jHcrnandcz et al. 2008) poses a serious problem for the 
paradigm that grain growth is a steady, monotonic pro- 
cess in disk evolution and planet formation. Addition- 
ally, small dust has been observed in the surface lay- 
ers of disks in clusters of different ages and environ- 
ments for hundreds of syste ms. The imp li cation s, as 
discussed mos t recen tly by lOliveira et al.l ()2010( ) and 
lOlofsson et al.l ()2010( ). is that small grains must be re- 
plenished by fragmentation of bigger grains, and that an 
equilibrium between grain growth and fragmentation is 
established. lOliveira et al.l ( 2010|) have shown that this 
equilibrium is maintained over a few million years, as 
long as the disks are optically thick, and is independent 
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of the population or environment studied. 

In this paper we present a comprehensive study of 
the mineralogical composition of disks around stars in 
young star-forming regions (where most stars are still 
surrounded by optically thick disks) and older clusters 
(where the majority of disks has already dissipated). 
Correlating the results on mean size and composition of 
dust grains per region, obtained in a homogeneous way 
using the same methodology, with the properties of small 
bodies in our own Solar System can put constraints on 
some of the processes responsible for disk evolution and 
planet formation. The Serpens Molecular Cloud, whose 
complete flux-limited Y SO population has b een observed 
by the IRS instrument (jOliveira et al.|[20lol ). is used as a 
prototype of a young star-forming region, together with 
Taurus, the best studied region to date. The sources 
that have retained their protoplanetary disks in the rj 
Chamaeleontis and Upper Scorpius clusters are used to 
probe the mineralogy in the older bin of disk evolution. 

Section [2] describes the YSO samples in the 4 regions 
mentioned. The Spitzer IRS observations and reduction 
are explained. The spectral decomposition method B2C 
(jpiofsson ct al. 2010) is briefly introduced in §[3l and its 
results for individual and mean cluster grain sizes and 
composition are shown in § 21 In § O the results are dis- 
cussed in the context of time evolution. There we demon- 
strate that no evolution is seen in either mean grain sizes 
or crystallinity fractions as clusters evolve from '-^l to 
8 Myr. The implications for disk formation and dissi- 
pation, and planet formation are discussed. In § [6] we 
present our conclusions. 

2. SPITZER IRS DATA 

The four regions presented here were chosen due to 
the availability of complete sets of IRS spectra of their 
IR-excess sources, while spanning a wide range of stellar 
characteristics, environment, mean ages and disk frac- 
tions (the disk fraction of Serpens is still unknown, see 
Table [J). 

The IRS spectra of a complete flux-limited sample 
of young stellar objects (YSO) i n the Serpens Molecu - 
lar Cloud have been presented bv lOliveira et al.l (|2010( ). 
based on program ID #30223 (PI: Pontoppidan) . As 
detailed there, the spectra were extracted from the ba- 
sic calibration data (BCD) using the reduction pipeline 
from the Spitzer Legacy Program "From M olecular Cores 
to Planet- Forming Disks" fc2d. iLahuis et a l. 2006i). A 
similarly large YSO sample ii i the Taurus star-fo rming 
region has been presented bv iFurlan et al.l ()2006l ). IRS 
spectra of all 18 members of the 77 Chamaeleon t is clus - 
ter were first shown by iSicilia-Aguilar et all ()2009[ ). 
while the spectra of 26 out of the 35 IR-excess sources 
in the Upp er Scorpius OB association were shown by 
iDahm fc Carpen ter (2009.) (the remaining 9 objects were 
not known at the time the observations were proposed). 
For the latter 3 regions, the post-BCD data were down- 
loaded from the SSC pipeline (version S18.4) and then 
extracted with the Spitzer IRS Custom Extraction soft- 
ware (SPICE, version 2.3) using the batch generic tem- 
plate for point sources. As a test, the IRS spectra of 
the YSOs in Serpens were also reduced using SPICE to 
ensure that both pipelines produce nearly identical re- 
sults. On visual inspection, no discrepancies were found 
between the results from the two pipelines, all objects 



showed the exact same features in both spectra. The 
similarity in outputs is such that the effects on the spec- 
tral decomposition results are within the cited error bars. 

Since the spectral decomposition method applied here 
aims to reproduce the silicate emission from dust parti- 
cles in circumstellar disks, the sample has been limited to 
spectra that show clear silicate features. The few sources 
with PAH emission have been excluded from the sample. 
PAH sources am ount to less than 8% in low-mass star - 
forming regions (|Geers et al.l[2006l : lOlI^ra et aLllMol) . 
Furthermore, spectra with very low signal-to-noise ratios 
(S/N) are excluded from the analysed sample in order 
to guarantee the quality of the results. In addition, for 
objects #114 and 137 in Serpens, and 04370-1-2559 and 
V955Tau in Taurus the warm component fit contributes 
to most of the spectrum, leaving very low fluxes to be 
fitted by the cold component. This produces large uncer- 
tainties in the cold component fit, and they are therefore 
not further used in the analysis. The low S/N objects 
rejected amount to less than 10 % of each of the Serpens 
and Taurus samples, so the statistical results derived here 
should not be affected by this removal. The final sam- 
ple of 139 sources analysed is composed of 60 objects in 
Serpens, 66 in Taurus, 9 objects in Upper Scorpius, and 
4 in 77 Chamaeleontis. The statistical uncertainties of 
the sp ectra were estimated as explained in lOlofsson et al.l 
((200a) . 

The great majority of the objects studied here are low- 
mass stars (spectral types K and M, see Table [5]). The 
study of mineralogical evolution across stellar mass is 
not the focus of this paper. Such a study would require 
a separate paper, in which the same techniques are used 
for low- and intermediate-mass stars. Thus, the statis- 
tical results derived in the following sections concern T 
Tauri stars, and not necessarily apply to intermediate- 
mass Herbig Ae/Be stars. 

3. SPECTRAL DECOMPOSITION AND THE B2C METHOD 

In order to reproduce the observed IRS spectra of these 
circumstellar disks the B2C decompositio n method, ex- 
plaine d in detail and tested extensively in Olofsson et al.l 
(|2010f l. is applied. Two dust grain populations, or com- 
ponents, at different temperatures (warm and cold) are 
used in the method, in addition to a continuum emission. 
The warm component reproduces the 10 /xm feature, 
while the cold component reproduces the non-negligible 
residuals at longer wavelengths, over the full spectral 
range (see Figure [T]). Each component, warm and cold, 
is the combination of five different dust species and three 
grain sizes for amorphous silicates or two grain sizes for 
crystalline silicates. 

The three amorphous species are silicates of olivine sto- 
ichiometry (MgFeSi04), silicates of pyroxene stoichiom- 
etry (MgFeSiOg), and sihca (Si02). The two crystalline 
species are both Mg-rich end members of the pyroxene 
and olivine groups, enstatite (MgSiOs) and forsterite 
(Mg2Si04). As further explained in lOlofsson et al.l 
(|201Q) . the theoretical opacities of the amorphous species 
are computed assuming homogeneous spheres (Mie the- 
ory) , while those for the crystalli ne species use th e distri- 
bution of hollow spheres (DHS, iMin et al]|2005[ ) theory 
so that irregularly shaped particles can be simulated. 

In addition, the three grain sizes used are 0.1, 1.5 and 
6.0 /zm, representing well the spectroscopic behaviour of 
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Fig. 1. — Example of the B2C modeling for object #15 in Ser- 
pens. The black line is the estimated continuum for this source. 
The red line is the fit to the warm component and the blue line is 
the fit to the cold component. The green line is the final fit to the 
entire spectrum. The original spectrum is shown in black with its 
uncertainties in light grey. 

very small, intermediate-sized and large grains. For the 
crystalline species, however, the code is limited to only 
2 grain sizes (0.1 and 1.5 /im). This restriction is im- 
posed because large crystalline grains are highly degen- 
erate with large amorphous gra ins (as can be seen in 
Figure 1 of Olofs son et al.l l201Clf ) , and because the pro- 
duction of large 6.0 /im pure crystals is not expected via 
thermal annealing (Gail 2004). 

The B2C method itself consists of three steps. First, 
the continuum is estimated and subtracted from the ob- 
served spectrum. The adopted continuum is built by 
using a power-law plus a black-body at temperature 
Tcont- The power-law represents the mid-IR tail of emis- 
sion from the star and inner disk rim. The black-body 
is designed to contribute at longer wavelengths, and is 
therefore constrained to be less than 150 K. Each dust 
component is then fitted separately to the continuum- 
subtracted spectrum. 

The second step is to fit the warm component to repro- 
duce the 10 /im silicate feature between '^7.5 and 13.5 
/im. This is done by summing up the 13 mass absorp- 
tion coefficients (A'spocios — 5, iVsizcs = 3 or 2, for amor- 
phous and crystalline species, respectively), multiplied 
by a black-body B^{Ty„) at a given warm temperature 
T 

The third step is to fit the residuals, mostly at longer 
wavelengths, over the entire spectral range (5 - 35 /im). 
This is done in a similar manner, for a given cold tem- 
perature Tc- The final fit is a sum of the three fits de- 
scribed, as can be seen in Figure [T] The entire fitting 
process is based on a Bayesian analysis, combined with 
a Monte Carlo Markov chain, in order to randomly ex- 
plore the space of free parameters. The resulting mean 
mass-average grain size is the sum of all sizes fitted, each 
size being weighted by their corresponding masses, as: 

where ai = 0.1 /im (small grains), 02 — 1.5 /im 
(intermediate-sized grains) and 03 = 6 /im (large grains) . 
Further details and tests of t he B2C procedure can be 
found in [olofsson et al.] ()2010D . That paper also demon- 
strates that the procedure is robust for statistical sam- 



ples, and that the relative comparisons between samples, 
which are the focus of this paper, should not suffer from 
the assumptions that enter in the procedure. The robust- 
ness of the procedure is evaluated by fitting synthetic 
spectra, and is discussed in detail in their Appendix A. 
The influence of the continuum estimate is also discussed, 
especially for the cold component for both grain sizes and 
crystallinity fractions, and it is shown that prescriptions 
that do not use large 6 /im grains (which are, to some 
degree, degenerate with the continuum) give flts that are 
not so good. 

For the amorphous grains, the B2C procedure uses the 
Mie scattering t heory to compute mass absorption coeffi- 
cients. However. iMin et all (|2007D found that they could 
best reproduce the extinction profile toward the galactic 
center using the DHS scattering theory, with a maximum 
filling factor of 0.7. The most striking difference between 
Mie and DHS mass absorption coefficients is seen for the 
0-Si-O bending mode around 20 /im. Here we investi- 
gate the influence of the use of DHS instead of Mie for 
amorphous grain with an olivine or pyroxene stoichiom- 
etry. We conducted tests on a sub-sample of 30 objects 
(15 in Serpens and 15 in Taurus). The conclusion of such 
tests is that it has a small influence on the quantities we 
discuss in this study. For the warm component of the 30 
objects, we find a change in the mean crystallinity frac- 
tion of -1.6% (the mean crystallinity for this sub-sample 
using DHS is 9.3% versus 10.9% using Mie), which is in 
the range of uncertainties claimed in this study. We also 
computed the mean slope of grain size distributions to 
gauge the effect of using DHS on grain sizes. On average, 
the grain size distribution indices are steeper by ~0.2 
(with a mean slope of -3.01 for this sub-sample using DHS 
versus -2.80 using Mie). Therefore, our main conclusions 
are preserved for the warm component. Concerning the 
cold component, the inferred crystallinity fraction using 
DHS is 22.5% versus 15.1% with Mie, a mean increase 
of 7.4%. For the mean slope of grain size distributions, 
a negligible decrease is found (-3.07 using DHS versus 
-3.01 for Mie). Again, the differences found are within 
our significant errors for the cold component and do not 
change any of our conclusions. 

It is important to note that the S/N generally degrades 
at longer wavelengths when compared to shorter wave- 
lengths. The lower S/N reflect on the cold component 
fits and will most likely result in larger uncertainties. We 
evaluate that the flts to the cold component are reliable 
and add important information on the dust mineralogy 
(albeit with larger uncertainties) and thus those results 
are included in the following discussion. 

4. RESULTS 

The IRS spectra of the 139 YSOs with IR excess dis- 
cussed in §[2] were fitted with the B2C spectral decompo- 
sition procedure. The relative abundances derived for all 
objects are shown in Appendix \X\ The S/N drops con- 
siderably for the long wavelength module of some of the 
objects studied (including all objects in Upper Scorpius 
and 77 Chamaeleontis) . For this reason, the cold com- 
ponent could not be satisfactorily fitted and no results 
for this component are presented for these sources (see 
Appendix . 

Due to the large number of objects, these results allow 
statistical studies on both the mineralogy and size dis- 
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tribution of the grains that compose the optically thin 
surface layers of disks in each cluster studied. The mean 
abundances of each species per region are presented in 
Table [21 where it can be seen that the majority of the 
dust studied is of amorphous form. In Table [3] the mean 
mass-average grain sizes and crystallinity fractions per 
region are shown. Mean sizes are in the range 1-3 ^m, 
without significant difference between regions. These re- 
sults are discussed in detail in the following sections. 

It is important to note that the comparison of results 
derived here for the different regions is valid because the 
same method, with exact same species, is used for all 
sources. The comparison of samples analyzed in distinct 
ways can lead to differences in results that do not cor- 
respond to real differences in composition. Nevertheless, 
in § 14.51 the results presented here are compared to liter- 
ature results for the same objects, when available, with 
generally good agreement. 

4.1. Grain Sizes 

The mean mass-averaged grain sizes for the warm 
((tiwarm)) and cold ((flcoid)) components are shown in 
Figure [2J for Serpens and Taurus (for the objects in Up- 
per Sco and ij Cha no results for the cold component are 
not available, see Appendix |A|. It is seen that the two 
clouds overlap greatly, and that the grain sizes derived 
from the different temperature components do not seem 
to correlate. To quantify this correlation, a Kendall r 
correlation coefficient can be computed together with its 
associated probability P (between and 1). r = 1(-1) 
defines a perfect correlation (anti-correlation) , and t — 
means that the datasets are completely independent. A 
small P, on the other hand, testifies to how tight the cor- 
relation is. For the warm and cold mean mass-averaged 
grain sizes for both clouds, r is found to be 0.14, with 
P = 0.07. This lack of correlation indicates that differ- 
ent processes are likely responsible for regulating the size 
distribution at different radii (iOlofsson et al...2010i) . 
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Fig. 2. — Mass-averaged mean grain sizes for the warm ((awarm)) 
and cold ((acold)) components. Black dots are the objects in Ser- 
pens, and red squares are the objects in Taurus. 

Although the average grain size in the warm compo- 
nent is bigger than that in the cold component within a 
given star-forming region, as shown in Table [H this dif- 
ference is mostly not significant. However, Figures [2] and 
[3] clearly show a difference between the range of grain 
sizes spanned in both components, with (acoid) never 
reaching near the biggest grain size modeled (6.0 fim) 
for any object. A possible explanation for larger grains 



at smaller radii, suggested by iSargent et al.l (|2009D . is 
that grains coagulate faster in the inner disk where dy- 
na mical timescale s are shorter. However as discussed 
by lOliveira et al.l (|2010( ) and in § 15.11 the mean dust 
size at the disk surface is not regulated by grain growth 
alone, but also by fragmentation and vertical mix. This 
means that faster coagulation at smaller radii cannot be 
uniquely responsible for bigger grains in the inner disk. 
Future modeling should try to understand this difference 
in mean grain sizes observed. 
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Fig. 3. — Distribution of mass-averaged mean grain sizes for the 
warm ((awarm), left panel) and cold ((flcold), right panel) compo- 
nents. r)ue to the low number statistics, the objects in Upper Sco 
and r] Cha have been merged together as an older cluster. 



Furthermore, Serpens and Taurus occupy an indistin- 
guishable locus in Figure [21 explicitly seen in Figure [3l 
A two sample Kolmogorov-Smirnov test (KS-test) was 
performed and the results show that the null hypothesis 
that the two distributions come from the same parent 
population cannot be rejected to any significance (14%). 
The older regions, although lacking statistical signifi- 
cance, show a distribution of mass-average grain sizes 
in the same range probed by the young star-forming re- 
gions (Figure [3]). This supports the evidence that the 
size distribution of the dust in the surface layers of disks 
is statis tically the same ind ependent of the population 
studied (|01iveira et al.ll2010[ ). 

The results here confirm those from IOlofsson et al.l 
([2010) that the mean differential grain size distributions 
slope for the three grain sizes considered are shallower 
than the reference MRN differential size distribution {a 
= -3.5). The mean grain size distributions slopes (a) for 
each region can be found in Table [3l 

4.2. Disk Geometry 

The amount of IR-cx cess in a disk is directly related 
to its geometry ( Kcnvo n fc Hartmannlll987l : iMeeus et al.l 
120011 : IPullemond et al.li200H ). Specificallv using the IRS 
spectra, disk geometry can be inferre d from the flux ra- 
tio between 30 and 13 /xm (F^n/Fi-i. prown et al.l 120071: 
lOliveira et al.ll2010t [Men'n et al.|[20ioD . A flared geome- 



try (1-5 < F30/F13 < 5), with considerable IR excess and 
small dust, allows the uppermost dust layers to intercept 
stellar light at both the inner and outer disk. For flat 



6 



Oliveira et al. 



disks < 1.5) with little IR excess, only the in- 

ner disk can easily intercept the stellar radiation as the 
outer disk is shadowed. Moreover, cold or transitional 
disks are interesting objects that present inner dust gaps 
or holes, producing a region with little or no near-IR ex- 
cess (5 < < 15). It is interesting to explore the 
effect of disk geometry on both the mean mass-average 
grain sizes and crystallinity fractions of the disks studied. 
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Fig. 4. — Top: Flaring index -Fso/i^is, used as a proxy for disk 
geometry, versus warm (left) and cold (right) mass-averaged mean 
grain sizes. Bottom: F30/F13 versus warm (left) and cold (right) 
crystallinity fractions. The YSOs in Serpens (black dots), Taurus 
(red squares). Upper Sco (blue stars), and 17 Cha (green triangles) 
are compared. 



not occurring at the same rate in the innermost regions as 
fu rther out in the disk. This is opposite to the conclusion 
of lWatson et al.l (|2009f) , who derive a correlation between 
inner and outer disk crystallinity from the simultaneous 
presence of the 11.3 and 33 /im features. The opacities 
of the crystalline species are more complex than those 
two features alon e, making the analys is here more com- 
plete than that of lWatson et al.l ()2009f ) . Our finding that 
the fraction of crystalline material in disk surfaces varies 
with radius can constrain some of the mechanisms for 
formation and distribution of crystals. 
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Fig. 5. — Crystalline fraction of the warm and cold components 
in Serpens (black dots) and Taurus (red squares). 



Figure [4] shows F30 / F13 as a proxy for disk geometry 
compared with the mean mass-averaged grain sizes and 
crystallinity fractions for both components and all re- 
gions studied here. No preferential grain size (correlation 
coefficient r = -0.14, P ^ 0.02, and r = 0.07, P = 0.33 
for warm and cold components, respectively) nor crys- 
tallinity fraction (t = 0.09, P = 0.10 for the warm, and 
r = -0.19, P = 0.01 for the cold component) is apparent 
for any given disk geometry. Similar scatter plots result 
for the mean mass-average grains sizes for only amor- 
phous (r = -0.12, P = 0.08 for the warm, and t = 0.13, 
P = 0.11 for the cold component), or only crystalline 
grains (r = 0.08, P — 0.17 for the warm, and r — -0.13, 
P = 0.10 for the cold component). Furthermore, no clear 
separation is seen between the different regions studied. 
The statistically relevant samples in Serpens and Tau- 
rus define a locus where the majority of the objects is 
located in each plot, which is followed by the lower num- 
ber statistics for older regions. Figure S] therefore shows 
not only that grain size and crystallinity fraction are not 
a function of disk geometry, but also that younger and 
older regions show similar distributions of those two pa- 
rameters. 

4.3. Crystallinity Fraction 

The crystallinity fractions derived from the warm and 
cold components (Cwarm and Ccoid, respectively) for Ser- 
pens and Taurus are show in Figure [S] No strong trend of 
warm and cold crystallinity fractions increasing together 
is seen (r = 0.10, with P = 0.10 for the entire sample). 
This fact implies that, if an unique process is responsible 
for the crystallization of dust at all radii, this process is 
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Fig. 6. — Distribution of crystalline fractions for Serpens (top), 
Taurus (middle), and Upper Sco and rj Cha combined (bottom). 
Similar distributions and the same range of fractions are seen for 
all clusters. 

A wider spread in crystallinity fraction is observed for 
the cold component than for the warm component (Fig- 
ure [6]), which is reflected in the mean crystallinity frac- 
tions for each sample (Table [3]). This discrepancy could 
be real, or an artifact due to the signal-to-noise ratio 
(S/N) being frequently lower at longer wavelengths (cold 
component) than that at shorter wavelengths (warm 
component), introducing a larger scatter. The differ- 
ence in Serpens is more significant ((Cwarm) — 11-0% 
and (Ccoid) - 17.5%). The left panel of Figure [7] shows 
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the cumulative fractions as functions of crystallinity frac- 
tions. Despite small differences between the warm (red 
line) and cold (blue line) components, the two cumula- 
tive fractions have similar behavior. If this difference 
is true, there is a small fraction of T Tauri disks with 
a higher cold (outer) than warm (inner) crystallinity 
fraction. This finding contrasts with that derived by 
Ivan Boekel et al.l ()2004D for the disks around 3 Herbig 
stars. Their spatially resolved observations infer higher 
crystallinity fractions in the inner than in the outer disks, 
albeit based on only 10 /im data. A larger sample of ob- 
jects with good S/N including both 10 and 20 /im data is 
needed to better constrain this point. In addition. Fig- 
ure |6] shows that younger and older clusters have similar 
distributions of crystallinity fractions. 
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Fig. 7. — Left: Cumulative fractions of the crystallinity fractions, 
for Serpens (dashed line) and Taurus (dotted line). Right: Cumu- 
lative fraction of the ration between the forsterite and enstatite 
fractions, for Serpens (dashed line) and Taurus (dotted line). The 
warm component is shown in red while the cold component is blue. 

As discussed by many authors, both the grain size 
and the degree of crystallinity affect the silicate features, 
therefore it is interesting to search for trends between 
these two parameters. In Figure [51 the mass-average 
grain sizes are compared to the crystallinity fraction for 
both warm (left panel) and cold (right panel) compo- 
nents. No obvious trends are seen in either component, 
neither any separation between region s . Thi s result sup- 
ports the discussion of lOlofsson et al.l (j2010D that what- 
ever processes govern the mean grain size and the crys- 
tallinity in disks, they are independent from each other. 
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Fig. 8. — Mass-averaged mean grain sizes versus the crystalline 
fraction for Serpens (black dots), Taurus (red squares), Upper Sco 
(blue stars), and rj Cha (green triangles). 



4.3.1. Enstatite vs. Forsterite 

The disk models of Gail (2004) consider chemical equi- 
librium of a mixture of solid and gas at high tempera- 
tures, allowing radial mixing of material. These mod- 
els predict a predominance of forsterite in the inner- 
most regions of the disk, while enstatite dominates at 



lower temperatures (being converted from forsterite). 
From the observat ional point of view, d ata on disks 
around T T auri (iB ouwman et aL| I2008D and Herbig 
Ae/Be stars (|Juhasz ct al. 2010) have shown the oppo- 
site trend: enstatite is more concentrated in the inner 
disk, while forsterite dominat es the colder, outer disk re- 
gion. IBouwman et al.l ()2008D interpret this result as a 
radial dependence of the species formation mechanisms, 
or a non-equilibrium of the conditions under which the 
species formed, contrary to the models assumptions. 

For the regions presented in this study, it can be seen in 
Table[5]for mean cluster values and in Tablc[5]for individ- 
ual objects that the results derived from thi s study gen- 
erally follow those of IBouwman et aTl (|2008D . with more 
enstatite in the warm component and, to a lesser extent, 
more forsterite in the cold component. The right panel 
of Figure [7| illustrates this for the cumulative fraction of 
the forsterite over enstatite ratios for individuals disks. 
However, this trend is not very significant given the un- 
certainties. 

4.4. The Silicate Strength- Shape Relation 

A correlation between the shape and the strength 
of the 10 /im silicate feature from di sks has been dis- 
cusse d extensively in the litera t ure (Ivan Boekel et all 
2001 iKessler-Si lacci e tahl l200l lOlofsson et al.l 120091: 
Pascucci et al. [1^009; Ol iveira et al. I I2010I: lOlofsson et al.l 
2010D. Synthetic 10 /im features generated for differ- 
ent grain sizes and compositions have been shown to fit 
well with observations, yielding grain size as the impor- 
tant parameter responsible for such a relationship. The 
degree of crystallinity of the dust also plays a role on 
the shape of this feature. However, as clearly shown for 
EX Lup (I Abraham et al.l 120091), and suppo rted by mod- 
els ()Min et al.l 120081: lOlofsson et al.ll2Ci09[ ). an increase 
in crystallinity fraction does not change the strength of 
the feature, even though its shape does change. Crys- 
tallinity is then understood as responsible for the scatter 
in the strength-shape relationship, and not the relation- 
ship itself. As a result, the strength and shape of the 10 
silicate feature yield the typical size of the grains in 
the upper layers of the dis k at a few AU from the star 
(|Kessler-Silacci et al.l [20071 ). The top panel of Figure [9] 
shows the results for Serpens, Taurus, Upper Sco and r] 
Cha. The bottom panel presents the median values per 
region, indicating the 15 - 85 percentil e ranges of the 
distributions. Overlaid are the models of lOlofsson et al.l 
((2009,) for different grain sizes (0.1 - 6.0 nm) gener- 
ated for amorphous silicates of olivine and pyroxene sto- 
ichiometry, and a 50:50 mixture. The difference in mean 
ages does not correspond to a significant difference in 
mean grain sizes between the different regions. 

With the mean grain sizes derived from the spectral 
decomposition, it is possible to further explore the valid- 
ity of using the strength of the 10 /im silicate feature to 
trace the sizes of grains in the surfaces of disks. The left 
panel of Figure [TUI shows the correlation between (awarm) 
and Sp^'^^ for all 4 samples. The Kendall t coefficient 
of -0.29, P — 0.01 supports the effectiveness of S'pea™ 
as a proxy for grain sizes, with smaller values of iSp^^™ 
implying larger grain sizes. 

On the other hand, it is also possible to test how the 
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Fig. 9. — Top: The ratio of normalized fluxes at 11.3 to 9.8 
lira (Sii.a/Sg.g) is plotted against the peak at 10 (•S'p^^^'™) 

for Serpens (black dots), Taurus (red squares), Upper Scorpius 
(blue stars), r] Chamaeleontis (green triangles). Bottom: Squares 
show the median values and crosses indicate the 15-85 percentile 
ranges of the distributions (top panel). Colored cur ves are derived 
from t heoretical opacities for different mixtures bv lOlofsson et ahl 
H2009l ^. The open circles correspond to different grain sizes, from 
left to right 6.25, 5.2, 4.3, 3.25, 2.7, 2.0, 1.5, 1.25, 1.0 and 0.1 /.tm. 




Fig. 10. — Left panel: Strength of the 10 fim silicate feature 
("^poak") versus the mass-averaged mean grain size for the warm 
component. Right panel: Strength of the 10 fim silicate feature 
versus crystalline fraction for the warm component. The best fit 
relationships are shown for reference. 

degree of crystallinity can influence the strength of the 
10 /im sihcate feature. The lack of correlation between 
Cwarm and S^^^^' (t = -0.07, P = 0.14), shown in the 
right panel of Figure [TU] for all samples, supports that 
the degree of crystallinity is not the dominant parameter 
setting the strength of the 10 /im s ilicate feature. These 
results argue against the results of iSargent et alj ()2009f ) 
that find a high crystallinity fraction and small grains 
fitting low strengths of the 10 /im silicate feature. Al- 
though it may be possible to fit a few spectra with a cer- 
tain prescription, a good model should be able to explain 
the robust relationship between the strength and shape 
of the 10 /xm silicate feature observed for large numbers 
of disks. Despite the many processes able to change the 
shape or the strength of this feature, only grain size has 
so far demonstrated capability to explain the observed 
trend. Our conclusion is that S^^'^^ and dust sizes are 
appropriately correlated. 



4.5. Comparison with other studies 

Dust composition results are available in the litera- 
ture for the disks in Taurus and 77 Cha (see Table H] for 
an overview). Sicilia-Aguilar et al.l ((2009il present their 
analysis in i] Cha considering the same 5 dust species 
and three grain sizes (enstatite in their model is the only 
species for which only the 2 smaller grain sizes are con- 
sidered), but for a distribution of temperatures derived 
using the Two Lay er Temperature Distribution (TLTD, 
lJuhasz et al.l I2OO90 decomposition procedure. For the 
same 4 objects, their mean amorphous fraction is 80.1 ± 
9.3 %. This result is consistent with the 82.8 ± 12.9 % 
mean amorphous fraction found here. The mean crys- 
talline fractions derived are 18.4 ± 10.7 % with TLTD 
a nd 17.1 ± 12.8 % der ived here. 

ISargent et al.1 (|2009D present their decomposition pro- 
cedure for 65 YSOs in Taurus. This method also takes 
into consideration a warm and a cold temperature, and 
makes use of two amorphous species (olivine and pyrox- 
ene) with two grain sizes (small and large), and 3 crys- 
talline species (enstatite, forsterite and crystalline silica) 
of a single size. Their mean warm amorphous fraction is 
82. 9± 19.3% and warm crystalline fraction is 17.1±19.3 
%, while here the derived fractions are 89.0 ± 6.6 % 
and 10.9 ± 6.6 % for the warm amorphous and crys- 
talline fractions, res pectively. For the cold component, 
ISargent et all ((2009) derive a mean cold amorphous frac- 
tion of 77.3 ± 19.9 % and cold crystalline fraction of 
22.6 ± 19.9 %, while here the values are 85.9 ± 10.6 
% and 13.9 ± 10.5 %, respectively. The consistently 
lower amorph o us (h igher crystalline) fractions found by 
ISargent et al.l ([2009) could he a result of their choice to 
use silica in crystalline rather than amorphous form (as 
used here). 

5. DISCUSSION 

5.1. Dust Characteristics 

Section |4] has shown that the disk populations in the 
four regions presented here, young and older, have very 
similar distributions in the two main dust parameters: 
grain size and composition. The large number of ob- 
jects in the two young regions studied occupy a region 
in parameter space of either grain size or crystallinity 
fraction that is also populated by the small number of 
older disks. The grain sizes derived for the cold compo- 
nent never reach the biggest grain size modeled (6 /xm), 
different from the warm component results that span the 
entire range in sizes. The crystallinity fraction does not 
seem to be correlated with mean grain size, warm or 
cold. Whatever processes are responsible for the crys- 
tallization of the initially amorphous grains, they should 
not only be independent from the processes that govern 
the grain size distribution, but they should also be able 
to work on bigger amorphous grains. Alternatively, the 
crystalline lattice should be able to keep itself regular 
during the coagulation of small crystalline dust to cre- 
ate big crystalline grains. The correlation between the 
strength of the 10 /xm feature and the mean grain size 
in disk surfaces, combined with the lack of correlation 
between crystallinity fraction and S^^^^, supports the 

wide usage of S^^l^^ a s a proxy for dust size in litera- 
ture (|van Boekel et al.ll2003t iKessler-Silacci et al.l 120061: 
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iPascucci et"aIll2009D. 

iBouwman et alj ()2008[ ) found a strong correlation be- 
tween disk geometry and the strength of the 10 /im sil- 
icate feature for a very small sample of T Tauri stars 
(7 disks), which points to flatter disks having shallower 
10 fim. features (i.e., big grains in the disk surface). 
Using results f rom s imila r decomposit i on pr ocedures, 
lOlofsson et all (|2010D and lJuhasz et all (|2010[ ) confirm 
this trend for larger samples of T Tauri (58 disks) and 
Herbig Ae/Be stars (45 disks), respectiv ely. Those trends 
are m uch weaker than that found by IBouwman et al.l 
(|2008f) . showing a larger spread. For the current even 
larger sample (139 disks), no significant trend is seen, 
indicating that the earlier small sample trends may have 
been affected by a few outliers . This result is similar 
to that found by lOliveira et all (|2010f ) for a large YSO 
sample 200 objects) using the strength of the 10 /^m 
silicate feature as a proxy for grain size (Figure 14 in 
that paper). As disc ussed by Oliveira et al ., the sedi- 
mentation models of iDullemond fc Dominilj ()2008() ex- 
pect a strong correlation of larger grains in flatter disks 
that is not seen. This means that sedimentation alone 
cannot be responsible for the distribution of mean grain 
sizes in the upper layers of protoplanetary disks around T 
Tauri stars. Furthermore, the lack of correlation between 
crystallinity frac tion and disk geomet ry is not in support 
of the results of I Watson et al.l ()2009l ) and iSargent et al.l 
(|2009t ). who find a link between increasing crystallinity 
fraction and dust sedimentation. 
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Fig. 11. — Left: Mean mass-average grain sizes vs. disk fraction. 
Serpens is not included because its disk fraction is not yet known. 
Filled circles represent mean warm grain sizes, and open triangles 
represent mean cold grain sizes. Error bars for the mean mass- 
average grain sizes are estimated using a Monte Carlo approach, 
sampling the errors of the individual objects. Right: Mean grain 
sizes vs. mean cluster age. Filled circles represent results for the 
warm component, while open triangles represent the cold compo- 
nent. The black points are YSOs in Serpens, red in Taurus, blue 
in Upper Sco and green in rj Cha. 



As discussed in lOliveira et al.l ()2010D for Serpens and 
Taurus, and confirmed by the addition of considerably 
older samples, there is no clear difference in the mean 
grain sizes in the disk surfaces with mean cluster age, 
which can be seen in Figure ill! This evidence sup- 
ports the discussion in that paper that the dust pop- 
ulation observed in the disk surface cannot be a result 
of a progressive, monotonic change of state from small 
amorphous grains, to large, more crystalline grains, or 
'grain growth and processing'. The fact that the distri- 
bution of grain sizes in the upper layers of disks does 
not change with cluster age implies that an equilib- 
rium of the processes of dust growth and fragmentation 
must exist, which also supports the existence of small 
grains in disks that are millions of years old whereas 



dust growth is a rapid proc ess ()Weidenschilling| 119801 : 
IDullemond fc Domini^ I2005D . That small dust is still 
seen in disks in older regions like Upper Sco and rj Cha 
argues that this equilibrium of processes is maintained 
for millions of years, as long as the disks are optically 
thick, but independent of them having a fiared or flatter 
geometry. 

5.2. Evolution of Crystallinity with Time? 

Literature studies of disk fractions of different YSO 
clusters with different mean ages show a trend of de- 
creasing disk fraction, i.e. disks di ssipating with time , 
over some few millions of years ()Haisch et al.l I2OOII: 
[Hernandez et al.l [2008) . This decrease is clearly con- 
firmed by the lower fraction of disks still present in the 
older regions studied here (Upper Sco and 77 Cha). Ac- 
cording to current planet formation theories, if giant 
planets are to be formed from gas rich disks, the op- 
tically thin, gas-poor disks in those older regions should 
already harbor (proto-)planets. Considering the evidence 
from small bodies in our own Solar System that suggest 
cons iderably higher cryst allinity fractions than IS M dust 
(see IWooden et all I2007i and .Pontoppidan fc BrearlevI 
l20T0l for reviews of latest results) , a crystallinity increase 
must occur. 
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Fig. 12. — Left: Crystallinity fraction vs. disk fraction. Serpens 
is not included because its disk fraction is not yet known. Filled 
circles represent mean warm crystallinity, and open triangles repre- 
sent mean cold crystallinity. Uncertainty for crystallinity fractions 
are estimated using a Monte Carlo approach, sampling the errors 
of the individual objects. Right: Crystallinity fraction vs. mean 
cluster age. Filled circles represent mean warm crystallinity, and 
open triangles represent mean cold crystallinity. The black points 
are YSOs in Serpens, red in Taurus, blue in Upper Sco and green 
in rj Cha. 



In Figure IT^ the mean crystallinity fraction per region 
is plotted against two evolutionary parameters: disk frac- 
tion (left) and mean age (right). Within the spread in 
individual fractions it is seen that, just as for grain sizes, 
there is no strong evidence of an increase of crystallinity 
fraction with either evolutionary parameter. This implies 
that there is no evolution in grain sizes or crystallinity 
fraction for the dust in the surface of disks over cluster 
ages in the range 1-8 Myr, as probed by the observa- 
tions presented here. Essentially, there is no change in 
these two parameters until the disks disperse. Starting 
from the assumption that initially the dust in protoplan- 
etary disks is of ISM origin (sub-/im in size and almost 
completely amorphous), it appears that a modest level 
of crystallinity is established in the disk surface early in 
the evolution (< 1 Myr) and then reaches some sort of 
steady state, irrespective of what is taking place in the 
disk midplane. Thus, the dust in the upper layers of 
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disks does not seem to be a good tracer of the evolution 
that is taking place in the disk interior, where dust is 
growing further for the formation of planetesimals and 
planets, at many times higher crystallinity fractions, to 
be consistent with evidence from Solar system bodies. 

If this is the case, within 1 Myr this surface dust must 
be crystallized to the observed fraction (~10 - 20 %). 
This result puts constraints on the formation of circum- 
stellar disks. One possibility is that this crystallization 
of the dust in disks mostly occurs during the embed- 
ded phase. In this early stage of star formation, where 
large quantities of material are still accreting towards 
the protostar, a fraction of the infalling material comes 
very close to the protostar and is heated to temperatures 
>800 K before it moves outwards in the disk. Alterna- 
tively, accretion shocks or episodic heating events could 
be responsible for thermally annealing the dust in the 
disk surface. 

The 2-D models of iVisser fc Dullemondl (|20Tot ) treat 
the radial evolution of crystals in time. According to 
these models, 100 % of the dust in the inner disk (< 1 
AU) is crystallized within 1 Myr. With time, the inner 
disk crystalline fraction drops as the disk spreads, and 
crystalline material is transported to outer parts of the 
disk. These models can help explain the rapid crystal- 
lization required to account for our results. However, the 
models do show a decrease in inner disk (< 1 AU) crys- 
tallinity fraction with time, which is not supported by our 
results. Since these models do not discriminate on ver- 
tical structure, but rather present crystallinity fractions 
that are integrated over all heights at a given radius, this 
decrease in crystallinity fraction is not necessarily con- 
nected to the surface of the disk. Thus the decrease in 
crystallinity fraction with time found in the models of 
IVisser & Dullcmond (2010) could be explained as a de- 
crease in crystallinity fraction just in the disk midplane 
where the bulk of the mass resides, but not in the surface 
layers, as our data indicate. That would imply that ra- 
dial mixing of these crystals is more efficient than vertical 
mixing, which is responsible for the crystallinity fraction 
decrease in the disk midplane.^ 

According to the models of iCieslal (|2007| ) for outward 
transport of high temperature materials, variations in 
radial transport dynamics with height produce vertical 
gradients in the crystalline fractions, such that the up- 
per layers of the disk will have lower crystallinity frac- 
tions than the midplane population. If that is the case, 
the observations discussed here, which probe the disk 
surface only, lead to lower limits on the real crystalline 
fraction of disk midplanes. In this scenario, planets 
(and comets) forming in the disk midplane would have 
higher crystalline abundances than those derived here 
for the disk surfaces, which are compatible with what 
has been observed in our Solar System. However, this 
model does not make predictions for the time evolution 
of the systems. Combining the vertical and radial mix- 
in g processes with evo l ution ary models such as those 
of IVisser fc Dullemondl (pfHOl) are needed to investigate 
whether older and younger disks could still show the 
same distribution of crystallinity fractions in the upper 
layers of disks, as observed here. 

6. CONCLUSIONS 



This paper presents the spectral decomposition of 
Sp itzer/IRS spect r a usin g the B2C decomposition model 
of lOlofsson et al.l ()2010D . Mineralogical compositions 
and size distributions of dust grains in the surface layers 
of protoplanetary disks are derived for 139 YSOs belong- 
ing to four young star clusters using the same method. 

Serpens and Taurus are used as prototypes of young 
regions, where most stars are still surrounded by disks, 
while Upper Sco and rj Cha represent the older bin of 
disk evolution, where a large fraction of the disks have 
already dissipated but some massive protoplanetary disks 
are left. The large number of objects analyzed allows 
statistical results that point to the main processes that 
affect the grain size distribution and composition of dust 
in protoplanetary disks. Furthermore, the usage of the 
same analysis method for regions of different mean ages 
allow a study of evolution of the dust parameters with 
time. 

Our large sample does not show a preferential grain 
size or crystallinity fraction with disk geometry, con- 
trary to earlier analyses based on smaller samples. Also, 
younger and older regions have very similar distribu- 
tions. The difference between mean mass-averaged grain 
sizes for the warm and cold components of a given star- 
forming region is small, however a considerable differ- 
ence is seen between the ranges of grain sizes spanned 
in both components. The cold mass-averaged grain sizes 
never reach the biggest size modelled (6 /im) while the 
warm mass-averaged grain sizes span the entire range 
of sizes modeled. The crystallinity fractions derived for 
inner (warm) and outer (cold) disks are typically 10 - 
20%, and not correlated. The cold crystallinity fraction 
shows a larger spread than the warm. No strong differ- 
ence is seen between the overall mean warm and cold 
crystallinity fraction. Within the crystalline dust popu- 
lation, more enstatite is found in the warm component 
and more forsterite in the cold component. The differ- 
ences are not very significant, however. 

The results of the spectral decomposition support the 
usage of the strength of the 10 /im silicate feature 
(■^pcakO as a proxy of the mean grain size of dust in the 
disk surface. This is supported by the correlation be- 
tween and mean grain size and lack of correlation 
with the mean crystallinity. 

Mean cluster ages and disk fractions are used as indi- 
cators of the evolutionary stage of the different popula- 
tions. Our results show that the different regions have 
similar distributions of mean grain sizes and crystallinity 
fractions regardless of the spread in mean ages of 1 - 8 
Myr. Thus, despite the fact that the majority of disks 
dissipate within a few Myr, the surface dust properties 
do not depend on age for those disks that have not yet 
dissipated in the 1-8 Myr range. This points to a rapid 
change in the composition and crystallinity of the dust in 
the early stages (< 1 Myr) that is maintained essentially 
until the disks dissipate. 

Astrochemistry at Leiden is supported by a Spinoza 
grant from the Netherlands Organization for Scien- 
tific Research (NWO) and by the Netherlands Research 
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TABLE 1 

Characteristics of the star-forming regions presented in this work. 



Region 


Dist (pc) 


Mean Age (Myr) 


Disk Fraction 


Ref. 


Serpens 


259-415 


2-6 




1, 2 


Taurus 


140 


~2 


~60% 


3, 4, 5 


Up Sco 


145 


~5 


-17% 


6, 7, 8 


T) Cha 


97 




~40% 


9, 10, 11 



Reference s: 1: Th e distance to Serpens is still under debate, different methods yield distances ranging from 259 fStraizys et al."1996') to 

415 pc IIDzib et al. 2010); 2: From Oliveira ct al. (200 9.), using d = 259 pc; 3 : From Kcnyon et al. (1994); 4: From Hartmann et al. 
(|200ll ): 5: From lTuhma n et al. (2010); 6: Fr om ide Zeeuw et al.l lfl999l): 7: From iBlaiTuwl 1119781): 8: FromiC arpentc r et al.l (12006|); 9: From 
iMamaick ct al. 1 119991) : 10: From lLuhman fc Steeghsl (120041 1; 11: From lMegeath et al . (2005). 



TABLE 2 

B2C mean composition of each star-forming region. 



Region 


Oli/Pyr" 


Silica 


Forsterite 


Enstatitc 




Warm Component 




Serpens 


81.3±11.7% 


7.8±6.5% 


5.8±4.9% 


5.2±4.9% 


Taurus 


79.4±9.4% 


9.6±7.0% 


4.4±3.3% 


6.5±4.5% 


Up Sco 


89.7±4.7% 


3.5±3.4% 


2.6±1.1% 


4.1±3.1% 


rj Cha 


75.0±13.4% 


7.8±4.6% 


6.8±4.0% 


10.3±9.8% 






Cold Component 




Serpens 


68.0±20.1% 


14.4±12.3% 


9.5±9.5% 


8.0±8.0% 


Taurus 


64.7±14.3% 


21.3±11.4% 


8.6±8.6% 


5.3±5.3% 



Amorphous olivine and pyroxene combined. 



TABLE 3 

B2C MEAN GRAIN SIZE AND CRYSTALLINITY PARAMETERS FOR EACH STAR-FORMING REGION. 



Region 


Number 


{^^warm ) 


(f^cold) 


(c^warm ) 


(acold) 


(Cwarm ) 


(C'cold) 






fim 


^m 






% 


% 


Serpens 


60 


2.9±1.3 


1.9±0.6 


-2.75±0.39 


-3.16±0.18 


11.0±6.9 


17.5±12.4 


Taurus 


66 


2.6±0.9 


2.4±0.6 


-2.83±0.31 


-3.02±0.15 


10.9±5.6 


13.9±10.1 


Up Sco 


9 


3.1±1.5 




-3.33±0.18 




6.8±3.3 




r] Cha 


4 


1.3±0.4 




-2.71±0.39 




17.1±10.6 





TABLE 4 

Comparison of mean mineralogical results from this study with literature studies. 



Region 



Taurus 
r] Cha 



This work Literature 
Amorphous Crystalline | Amorphous Crystalline 



89.0±6.6 % 
82.8±12.9 % 



Warm Component 



10.9±6.6 % 
17.1±12.8 % 



Taurus I 85.9±10.6 % 



82.9±19.3 % 17.1±19.3 % " 

80.1±9.3 % 18.4±10.7 % 
Cold Component 

13.9±10.5 % I 77.3±19.9 % 22.6±19.9 % " 



" ISareent erall (|2009|) ; '' ISicilia-Aeuilar et al.l ( l2009t) 
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APPENDIX 

RELATIVE ABUNDANCES OF SPECIES 

The relative abundances, as resulting from the B2C compositional fitting (§[3]) to the IRS spectra of protoplanetary 
disks in Serpens, Taurus, Upper Sco and r/ Cha are shown in Table [51 Since the opacities of amorphous olivine and 
pyroxene are degenerate, the abundances of these two species have been added into one, marked 'Oli/Pyr' in the Table. 
Furthermore, 'Sil' designates the amorphous silica, and the crystalline enstatite and forsterite are marked as 'Ens' and 
'For', respectively. In the table, the first line of a given object corresponds to the results of the fit to the warm 
component and the second line to the results of the cold component. For some objects (20 in Serpens, 28 in Taurus, 
all in Upper Sco and f] Cha), the S/N drops considerably at longer wavelengths and the results of the procedure are 
no longer reliable. For these sources, the cold component could not be fitted satisfactorily and, in Table [SJ only the 
warm component results are shown. 
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TABLE 5 

Dust composition derived using the "B2C" procedure.'' 



ii:> spT 


Oli T'm" '/ 


■ : 




'-■il '-5 


Oli ryi 


En;- 


1^1 .1 


'-■ii '-5 


( )ii Vx r '-5 


sil % 




^U,1 ,.m) 


(11,1 ,,m 


Ul.l ;/ni) 


(11,1 ,jmj 


{IJ, ,,ni 


(1,:, ,,m) 


(l,,--, (im) 


(1,,:, ,,m) 


((),(( ;<mj 


((,,() ,jm) 



Serpens 



1 


K2 


— 15.3 


1-0+nS 
— 0.4 


2.6 + 2-0 
— 1.4 


0-1+?-? 

— 0.0 


53.7+16^5 
— 23.0 


3 5+3-5 
' —2.4 


0-0+?-? 

— 0.0 


l-0+?-n 

— 0.0 


7 8+3-8 
' —0.0 


0.6+1-? 
— 0.0 






38, 4+1^-5 


0.0+1-0 

— 0.0 


8.7+5-5 

— 3.6 


1-2 + ?? 
— 0.0 


29.3 + 6,2 
— 15.3 


6.8+^3 
— 2.8 


0.2 + J^l 


0.3 + 1-1 
— 0.0 


9-7+11 


^'^—1.6 


3 


MO 


— 7.3 


0.0+2-f; 

— 0.0 


0-9+^i 

— 0.4 


0.0 + 1-7 
— 0.0 


43.9 + ;0^8 
— 18.6 


— 0.2 


2.7+3-8 


0-0+?? 
— 0.0 


15.7+3-2 
— 2.8 


0.0+10 

— 0.0 






— 8.6 


o.o+io„5 

— 0.0 


10.2+iV 

— 4.9 


0.0+16^3 
— 4.U 


32.9+34^1 


4.9+19,2 
— 5.5 


0.0+13-9 
— 2.9 


12.7+22/ 
— 5.6 


10.2+7-1 

— 4.7 


1-8+?-? 
— 0.6 


6 


K5 


25.4+^-^ 

— O.U 


0.0+0-'; 

— u.o 


0.8+1-1 
— U.O 


0.0+0-3 

— U.O 


38.2+11^6 


i-o+nn 

— 0.0 


lO+nn 

— 0.0 


0.0+1-1 

— 0.0 


28.8+1-5 
— O.U 


4.8+3-6 
— 1.8 






6.4+3-8 
— 3.3 


4.1 + 2-1 

— 2.3 


5-1+nn 
— 0.0 


0-5+;,;, 

— 0.0 


46.1 + 19-8 
— 2U.3 


4-7+?-? 
— O.U 


1-2+?? 
— U.O 


2-7+?? 

— 0.0 


25.4+8-6 
— 7.6 


3-7+?? 

— O.U 


7 


MO 


— 0.9 


0.0+2-9 

— 0.0 


— 0.0 


0-0+?? 
— 0.0 


75.6 + 26-0 

— 23.9 


0.0+5-6 
— 0.1 


3.7+*-7 

— 3.0 


0.0+5-1 
— 0.0 


0.0+9-8 

— 0.0 


le.i+t,!. 

— 13.7 






— 0.5 


7-2+6? 

— D.y 


8-1 + 6? 
— 6.0 


14.1 + 8-3 
— 11.5 


6.0+7-7 


14.1 + 8-8 

— 11.7 


0.0+3-0 

— 0.1 


24.7+1,1-3 


25,6+11-0 

— 13.9 


0.0+3-5 


9 




5g 3+10.9 


o.o+i-f, 

— 0.0 


— 0.0 


0.0+0-7 

— 0.0 


27.4+5-7 
— 7.2 


3.8+2-t 

— 2.1 


2-5+0? 

— 2.0 


0.0 + 1-3 
— 0.0 


3,6+1-1 

— 0-0 


+ 17 
^' ''' — 0.0 






24,5+1^3 


— 0.0 


14.4+13/ 
— 9.8 


16-9+?i? 
— 15.0 


9.4+6-7 


0-0+^1 
— 0.8 


0-0 + 00 


7.5 + 5-1 


24,4ll/;| 


2.7+7-2 

— 3.2 


10 




1 g + 5,9 


0-0+00 


2.3+2-1 


4.3 + 2^J 


37.6+17^3 


7.2+1-2 


5.3 + 6-6 


0-2 + ?-? 


31,6 + /,5 
— 8,3 


9.6+1-6 






37,4+f -1 
— 9.4 


2.0+10-1 
— 0.0 


0.0 + 0-'2 
— 0.0 


0.0+19-4 

— 0.0 


.39.0 + 2-3-'6 
— 10.3 


0.0 + 20-1 
— 0.0 


3.3 + 9-6 


0.0 + 9-3 
— 0.0 


14,8 + 3521 


3.6+11,1 
— 0.5 


14 


M2 


ii,i+*V 


0,0+0-2 


3.0 + 7-1 


0.0 + 6-9 


-54.7+-3J^3 


4.8+200" 


0-0+?-? 


0.0 + 8-7 
— 0.0 


13,5 + 0-3 
— 5.8 


12.8+// 

— 4.0 


15 




— 7.0 


0-2 + i^-^ 
— 0.0 


2.9 + 1-0 
— 1.3 


0.0+0-0 

— 0.0 


51-3 + ?«3, 
— 16.1 


3.0+3-t 
— 1.5 


0-0 + 1? 
— 0.0 


0.8 + 1-8 

— 0.0 


1-2+?-? 
— 0.0 


5-0+?-? 

— 0.5 






— 12.9 


2.4 + }-? 
— 1.6 


2-2 + ^-^ 
— 1.0 


1 1 + 2-8 
— 0.9 


34.9 + 9^9^ 


00+?? 

— 0.0 


5,4+2-5 


0.0 + 2-1 
— 0.6 


14,6+1-7 

— 4.6 


3.0+5-1 

— 1.9 


21 




2,2 + = -" 
— 0.0 


0,0+1-0 

— 0.0 


2-7+?-n 
— 0.0 


0-1+?? 

— 0.0 


20.2+;0-'l 
— 12.9 


01+n? 

— 0.0 


19,2 + 11,-1 
— 9.1 


0,0+1-2 

— 0,0 


43,3+:6,--5 
— 2 1 - 3 


12.3+!-? 

— 9.1 






33.0+111 


— 0.0 


0.0 + 2-? 
— 0.0 


1 3+11-0 
— 0.0 


42.5 + 22-6 


2-7+;,v 

— 0.0 


0-5 + ?-? 


0,0+^1/ 


9,6 + 8^7 


8.7+5-2 


29 


M2 


— O.U 


0.0+3-1 
— 0.0 


1-0+?;? 

— 0.0 


2.7+6-6 
— 0.3 


71.0 + 16-0 


0-0+?? 
— 0.0 


4,4 + 0-4 


0.0+7-9 
— 0.0 


13,0 + fl 


5.6+11/ 


30 


Ml 


— O.U 


0-1+no 

— 0.0 


3-4+^1 


4-5 + ?1 
— 0.2 


12.5 + 8-5 
— 6.4 


12.5+8^1 
— 3.4 


1 3 + 5-9 
— 0.9 


1 5 + 5-0 
— 0.8 


64,2ll8;J 


0.0+21/ 

— 0.0 






— 13.8 


1.7+ IV 
— 3.2 


5-5+lV 
— 4.3 


8.1+2I36 

— 6.3 


8.0 + 58-7 
— 19.4 


9.i + i,*;9 


29-6+"-? 
— 24. 3 


11.4+37-5 

— 10.3 


14,0+36-3 
— 10.3 


2.1+9-6 

— 2.3 


31 




1,2 + 0-8 
— 0.0 


0.0+0-8 

— 0.0 


0.8+0-5 
— 0.0 


0.0+1? 

— 0.0 


15.0 + 2-6 


-j^5. 1^ 
^■^ — 0.0 


3.4+15 
— 1.3 


0-1+?-? 

— 0.0 


76,1 + ^-3 

— 5.4 


06+?0 
— O.U 






21.8 + ^-^ 
— 4.9 


— 0.0 


0.0+0-7 
— 0.0 


16,7+8-7 


24-n+l-t 

— 4.0 


2.3 + 2^0 
— 0.0 


5+1^7 
— 0.0 


5.4+1-1 
— 2.6 


9,2+1-8 

— 5-2 


17.6+^-^ 
— 4.7 


36 


K5 


22. 3'''^-^ 
— 6.2 


0.9+1-6 
— 0.0 


— 0.0 


0.0+0^8 
— 0.0 


51.9+;2-l 
— 12.3 


i^o+?-? 

— 0.6 


3.3 + 2^? 
— 2.5 


0.0+1-6 

— 0.0 


10. 8 + ^1 
— 2.6 


9.0+5-7 

— 4.8 






37.4111-0 


0.0 + 3-0 
— 0.0 


1.6 + 2-7 
— 1.0 


2.9 + 8-6 


38.i+;i^o 


0.0 + 2^4 


0.0 + 2^9 


0.2+8-7 


19,6+5-0 


0,1+0-3 


40 


M7 


0,0 + 8-3 
— 0,0 


0,0 + 8-4 
— 0.0 


2-9+n8 
— 0.8 


10.9 + J-6 


41.2+17'3 
— 25. 


22-5+ll„0 
— 8.8 


17.9+11/ 
— 7.4 


0,0 + 5-5 
— 0.0 


4,5+9-7 
— 0-2 


0-0+?-'? 

— 0.0 






24,2+ 


16,6+llf 


3 j + 16.4 


0-0+n? 

— 0.0 


18.2+O35 

— 13.8 


2,1 + 12,8 


13,4 + 5^8 


0,0 + 3-0 
— 0.0 


0-2+?-? 

— 0-0 


22,3+51-2 


41 


K2 


— 0.0 


6,3 + 10„l 
— 0.0 


3.9 + 10-" 
— 0.0 


0-8 + ?? 
— 0.0 


^ + 6.8 

66.6 KO n 

— 58.2 


3.4+16-,3 
— 0.0 


16,4 + 11/ 
— 8.6 


0-0+?? 

— 0.0 


2,6 + // 
— 0.0 


0-0+nV 

— 0.0 


43 


MO. 5 


0.0+13„0 
— 0.0 


O-O+n'n 
— 0.0 


4-2+^7; 

— 0.0 


4-6+?-n 
— 0.0 


42.1 + 19-1 
— 27.0 


7.1 + 8,-9 
— 5.1 


2.7+7-1 

— 2.0 


0-0+?? 

— 0.0 


37,2+l?-3 
— 19.3 


2.0 + // 
— 0.0 


48 


M5.5 


7.0+1'-^= 


0-1+nn 

— 0.0 


2.0 + 10„2 
— 0.0 


0.6+?-? 
— 0.6 


73.3 + 28-2 

— 35.2 


0.0+9-2 

— 0.0 


0.0 + 13/ 
— 0.1 


1-1+?-? 

— 0.7 


6,0 + ?^t 

— 2-1 


9.9 + 1-6 

— 6-9 






46,0 + 23-0 

— 40.6 


+ 10.9 
' ■ —3.3 


2.5+10j5 


0.3+2I33 


18.7+;0-0 
— 17.6 


0.0 + n" 1 
— 0.1 


4.1 + 12/ 


8.6+5-8 

— 8.1 


71+4-? 


7 + 

—6-6 


53 


M2.5 


0.0+3^ 
— 0.0 


5.0+20^0 

— 0.0 


5.4 + 18^5 
— 0.0 


5.9+20^1 

— 0.0 


30.2+18-| 


10.1 + 10/ 


0.0+11/ 


2.3 + 15/ 
— 0.7 


32.9 + 27/ 
— 6.7 


s,i+o,-i 

— 5.3 


55 


K2 


16.9+11„3 
— 4.0 


l-O+n'n 
— 0.0 


2-3+?i 
— 0.9 


o.o+i^i 

— 0.0 


71.0+16-2 
— 26.2 


0.0 + 3-0 
— 0.0 


0.0 + 2-7 
— 0.0 


1-9+?? 

— 0.0 


12+?-? 

— 0.8 


2,1+0-? 

— 0.7 






27,7+81-3 


4.8+fiO 


21.5+34-7 


0.0+22„0 

— - 


20.5 + 7J-1 


15-7+72„9 
— 9.9 


0.0 + 22-3 


4.8 + 65/ 


2.9 + 5^? 


1,9+5-0 

— 1.8 


56 




7,8 + 8-3 


1.7+1-0 


2.5+2-4 


1. 1 + 3^4 


44.6+16-6 
— 28. 5 


5.0+7-9 


1-7+7-3 


0,0 + 1-6 


29.5 + ;i^8 


6,l + /,2 
— 0.1 






— 17.0 


1.6 + f'l 
— 1.0 


0.0+0-1 

— 0.0 


2.2 + 2 V 

— 0.0 


29.8+;i^? 
— 12.0 


0.0+15o8 
— 0.0 


1 1 + 7--5 
■ —0.0 


0-0+?-? 

— 0.0 


19.3 + 21/ 
— 7.9 


1.6 + ?7/ 
— 0.0 


57 




— 0.0 


0.0+0-1 

— 0.0 


5.7+7-2 


3-0 + 11 


72.0 + 29^0 
— 40.9 


10.5+13.3 
— 6.6 


0.0 + 7-1 
— 0.0 


0-0+?-? 

— 0.0 


6.5 + /.1 
— 0.0 


0.0+9-5 
— 0.0 






23.0+™„0 
— o.o 


— O.U 


14.4+tO-O 


0.7+131-2 
— 0.0 


13.1+55„6 
— 4.8 


— 0.0 


4.5+61/ 
— 0.0 


1.7+f,5 
— 4.5 


25.ll|5/ 


17.2 + 81/ 
— 5.0 


58 


K7 


6.6 n c 

— 0.5 


0.0+1-0 

— 0.0 


0.6+1-i 

— 0.0 


1-2+?;; 

— U.O 


29.9+^6 
— 6.7 


1-0+?? 
— 0.7 


7.1+1-0 


0.0+1? 

— 0.0 


28.3+6-5 
— 11.5 


25.4+6-1 
— 8.7 






32.2+16.8 
— 16.4 


3.2+15„6 
— 0.0 


0.0+3-7 
— 0.0 


— 0.0 


57.0+24-9 
— 36.3 


0.0+7-5 
— O.U 


0.0+8-0 
— U.O 


0.0+6-2 
— 0.0 


4.7+5-9 


2.8+2-5 

— O.U 


60 


MO. 5 


o.o+io„i 

— 0.0 


9-8+nV 

— 0.0 


5-2+;,^ 


0.0+3-9 

— 0.0 


41-7+IV7 

— 21 .7 


8.5+13,1 


3.3+1*/ 

— 2.3 


0.0+13.7 

— 0.0 


4.9+;i/ 

— 1.7 


26.5+57-0 

— 14.6 


61 


MO 


— 0,0 


0.0+0-0 

— 0.0 


+13 

"'""^ — 0.0 


0-0+nn 
— 0.0 


^^■■'^ — 13.2 


"+5 1 
— O.U 


6.2 + 3-5 

— 3.0 


o,o+?-i 

— 0.0 


11-3+3? 

— 3.0 


13.,3+1-? 
— 6.3 






3S,2 + |-^ 


0.0 + 0-8 


0-1 + 0? 
— 0.0 


0.6+1^7 
— 0.0 


38.7+6-6 

— 6.3 


0.0+1-6 
— 0.0 


1.1 + 09 


1,2+2-5 


7 5 + 5-7 

— 0-0 


12.7 + 1-5 


65 




0,0 + 2-^ 
— 0.0 


0.0+2-0 
— 0.0 


6.0+2-? 
— 2.5 


0.0+1-0 

— 0.0 


21.3 + 9^9 
— 16.6 


8-6+",-? 
— 5.4 


8.3 + ^-2 
— 3.2 


0.0+1-3 
— 2.7 


41.8+15-0 
— 22.4 


i3.9+^;p . 

— 10.4 






54,9 + ?0-? 
— 19-1 


0-1 + nn 
— 0.0 


0.0+0-3 

— 0.0 


8.8+16„2 

— 6-6 


6.2+9-7 

— 2.8 


0-3+?? 

— 0.0 


0.0+1,-0 

— 0.0 


2.7+3-8 

— 1.4 


8.5 + 28/ 

— 2.5 


18.4+20/ 
— 7.7 


66 


K5 


— 0.0 


0.0+1-0 

— 0.0 


0-1+nn 

— 0.0 


0-8+?-n 
— 0.0 


46.5 1 K f. 
— 15.6 


08+?n 
— 0.0 


2.9+5-? 

— 1.6 


0.0+2-3 

— 0.0 


38.2 + ;^-^ 
— 18-7 


8,7+!-2 

— 2.5 






33,5+;0-5 

— 16.5 


3.9+10„l 
— 0.0 


4-2+iV 
— 0.0 


0-0 + ?? 
— 0.0 


26.3+25„8 
— 8.9 


1^7+nV 
— 0.0 


0-4 + ?-? 
— 0.0 


3.2+// 

— 0.0 


ii.o+?i/ 

— 0.0 


15,7+6-8 


71 


M3 


6,4+28^2 
— 0.0 


0.0 + 1-8 
— 0.0 


2.3 + 0-8 
— 0.0 


0-8 + ?-? 
— 0.0 


24.2 + ?f-6 
— 11.4 


3.0+2031 

— 0.3 


4.7 + 2-5/ 
— 1.6 


1,9+12,1 
— 0.5 


56. 6+55-6 
— 10-9 


o,o+/oi 

— 0.0 


74 




0.0+nn 
— 0.0 


0.0+1-7 

— 0.0 


6.3+3-7 
— 2.8 


13-4+1? 
— 6.3 


0-0+?-? 

— U.O 


5-8+?1 
— 3.0 


16.1+1 

— 10.0 


0.0+1-7 

— 0.0 


58.3+15^0 


0,0 + // 
— 0.0 






50.7+^^-t 
— 19.5 


— 0.0 


0-l+?n 
— 0.0 


0-0 + ?V 

— 0-0 


24.2+26.5 
— 9.8 


1.8 + ?3-,4 
— 0.0 


5-4+?V 
— 0.0 


0.0+12/ 
— 0.0 


9.9+2,2/ 

— 5.9 


7.8+23,-l 
— 7.1 


75 




o.o+';2„5 

— 0,0 


0.0+10.2 
— 0.0 


4.1 + 28-5 


0.1+20^2 


15.0+>15-0 
— 0.0 


5.2 + 36o° 
— 0.0 


0.0+?2/ 
— 0.0 


o.o+?i/ 

— 0.0 


59.6+12^6 


16.0 + 69/ 
— 0.0 


76 


Ml 


5 3 + 0-'' 


0-0+00 


2.613-3 


1.9 + -3-1 


65-0 29 Q 


0-0+?-? 


1.3 + 5-3 


0-0+?? 


23,7l'l3-5 


01 + ?-? 

— 0.0 






41,8+1^5 
— 23.6 


3.6 + *-5 
— 1.9 


0.0 + ^-^ 
— 0.0 


0.0 + 2-8 
— 0.0 


29.6+;2-'4 
— 12.1 


0-0+?-? 

— 0.0 


12-9+"7-2 
— 5.2 


0.0 + 2-1 
— 0.0 


12.2 + ^^'-^ 
— 4.3 


0.0+5-1 

— 1.0 


80 




0.3+1-2 
— 0.0 


0.0+0-2 
—0.0 


1-2 + ?'? 
— 0.0 


0.0+0-1 


7 5+2-0 

— 0.0 


3.2+;-3 


3-1+1-0 


0,0 + 0-6 


84.7+8-? 

— 6.D 


0.0+1-; 






6 4+50-4 
— 0.0 


a ,+29.8 
^-•^-0 


6+12-* 


7-4l?V 


15 4 + 50.2 


6.2i;,8/ 


0.8l// 


3.1 + 0-0 


22.6178/ 


29.2167/ 


81 


M5 


10 7+^-^ 


+ 0-8 


+ 0-4 


4 3 + 1-2 


54.5 + 5.4 


+ 0-7 
"-"-0.0 


6 1 + 1-9 
"1-2.0 


0+0-5 
"-"-0.0 


15 6+2-5 


8-8i5-:i 








0+1-1 


0+0-* 


0-oto,? 


31.418-1 


7-0l^2-.*2 


+ 0-3 
"-"-0.0 


1-01^2-,? 


7.715-0 


l-5t?-.^ 


86 


M5.5 


i5.7t?:i 


0+0-5 
o-"-o.o 


8+0-7 
0-8-0.0 


0-ot?:? 


36-01^;^^ 


o-ol?,? 


3 3+1-1 


o-oi?:? 


ll-lll-.l 


1 + 0.7 
"1-0.0 






96 7 + 83. 1 


19 r)+ 28 . 5 

12--'-0.0 


0+38-4 
0-0-0.0 


1-oiJV 


18.2 + 55^5 


10 + 27-7 
l"-"-0,0 


5 + 20,3 

"-"-0,0 


6 + 29-6 
"-"-0,0 


10 7+39,6 
10- '-0,0 


,3 7+38.7 
1-'- '-0.0 


88 


MO. 5 


0+5-0 

0-0-0.0 


0+10 
o-"-o.o 


8 6+4-4 
8-0-0.0 


+ 0-4 
"-"-0,0 


67 1 + 20-5 
"'-1-17.3 


10 3+7-1 
l"-3-5,0 


+ 3-6 
"■"-0,0 


+ 2,4 
•'■1-2,0 


9 9+6-3 


9+1-2 
"--'-0-6 






6 + 30-9 
°-O-10.0 


7+37.8 
"■'-9.0 


33 3+22.2 
•^•^-^-13.7 


6 0+31-7 
6-"-9.2 


2.1 + 21,3 
— 6.5 


0+7-7 
"■"-0,0 


on + 29,1 
30-6li2.4 


6 + 0^9 
"■"-3,0 


17,6152.8 


3 1 + 11-2 
3-1-3.5 


89 


K5 


0+0-^ 

O-"-0.0 


0+0-2 

o-"-o.o 


2 6+10 


1 6 + 0-7 
1-6-0.0 


0+0-7 
"-"-0.0 


0+0^6 
"■"-0.0 


1 4+0-9 
1-1-0.0 


0+06 
"■"-0,0 


94.4l/,-8 


o-ol?-7o 


90 




26.0+20-0 


o+'i-i 

0-"-0.1 


3 + 4.4 
■'■"-1,6 


+ 2-7 
"■"-0.0 


49 1 + 24.7 
*9-i-30.7 


3 7+9-5 
3-'-1.8 


+ 10-1 

"-"-0.0 


2,8l-5;5 


K 4 + 16.4 
■'.1-1,2 


10.0152/ 






o, ,+65:4 

31-1-0.0 


7 5 + 13-4 


+ 51,:! 

8--'-(),0 


11 7^21.3 
"■ ' -0-0 


16 R+54.1 
16-8_o.o 


3 6 + 38-8 
3-6-0.0 


9 2+58-7 
9-2-0.0 


+ 24,8 
-(1,0 


12 5 + 52.2 

i^.''-o,o 


q 2 + 39,3 
-'■2-0,0 


92 


MO 


1 1 + 0-7 

1-1-0.0 


0+0-9 

o-"-o.o 


4 s + 4^8 
^■8-0,0 


2 3 + 5.3 


35-8_8 5 


5 7+7-5 
5-7-2.1 


+ 3-5 
"-"-0.0 


+ 2^1 
"■"-0,0 


44,0155.1 


6 3+11^1 
63-2,4 






19.6^0-7 


n s+3-4 
0-8-0.0 


2 + 0^0 
"■2-0.0 


+ 4,6 
■'■"-0,0 


31 2 + 0-3 
31-2-15.1 


4+6-5 
"■ -0.0 


6 8+3-8 
6-8-3.9 


+ 10^0 
"■"-0.0 


37,9i5o;i 


0+12^6 
"■"-0,0 


96 


Ml 


18-2i?^8 


9 c + 3.8 
2-'>-2.0 


2 8+''^6 
28-1. 6 


0+5^9 


56.5123-3 


3 8+12-9 
3-^-1.8 


0+15-2 
"-"-0.0 


3 0+9^2 
3-"_o.8 


7-i±l:l 


5 8 + 10-1 
5-8-1.5 



Amorphous olivine and pyroxene combined. 

For each object, the first line corresponds to the warm component abundances, while the second line corresponds to the cold component abundances, when available. 
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TABLE 6 
Continuation 





((J.l ,,1,1) 


(0.1 ,,111) 


(0.1 ,11,1) 


((J.l ,1111) 




(I..-1 ,1111) 


{i.r, ,,111) 


(i.ri ,1111) 


((,.() ,,111) 


(0.0 ,,111) 


100 






o+i^-" 

"•"-0.0 


2 s+*2-3 
2-»-0.0 


"-°-o.o 




q = + 81.9 

^•°-o.o 


+ 3*-l 
"■"-0.0 


, ,+lU6.2 
1^^-0.0 


8i.8l-;3 


3 , +11.6 
3-1-4.1 


101 




,+26.6 
''^■'--6.9 


o+^-i 

"-"-0.0 


2-"-1.4 


1 0+*'-2 
^■"-0.0 


60.2I232.6 


1 Q+8-8 


2 6+''-5 


0+53 
"■"-0.0 


6.6lJ0^1 


1 6 + 5-2 
1-6-0.0 








+ 2-7 
"■"-0.1 


2 + 2-7 
"-2-0.0 


2 5 + 3-1 
^-"-1.9 


55.9l^^lg 


0+3-7 
"-"-0.0 


9 5 + ''-8 
"-"-6.1 


0+2 7 
"■"-0.0 


13.6 + 3^5 


o-i^?:S 


103 




"•"-0.0 


0+0^3 
"•"-0.0 


2 4 + 9-6 
2-*-0.0 


3 + 9-9 
3-^^-0.0 


30 n+7-0 
3»-0-5.0 


2 0+2-3 
2-"-1.2 


7 8+2-4 
7-8-2.1 


0+18 
"■"-0.0 


27.6l6,0 


18-911:? 






11 0+^''' 


5 S + '' 2 


0+3-9 
"■"-0.0 


5 8 + 3-3 


3_ ,+18.1 
3^3-^-20.5 


0+6-3 
"-"-0.0 


5-ii^o:? 


^■2l"o:!; 


15.7l8'5 


17 3+2^3 
^'■3-17. 7 


104 




'^"-l.l 


9 + 61 
"•"-0.2 


3 7+f'^2 

3'-2.3 


4 3+*-l 
*-3-1.8 


70 4+35.4 
'"■*-48.S 






"■olSj.i 


14.2117^3 


1 + 14-8 
"■^-0.0 


106 


M3 


0+16 
"■"-0.0 


+ 23.7 
•^•1-0.0 


4 n + 17^1 
^■"-0.0 


= -71^0^0= 


1 + 24-9 
"•^-0.0 


10.012,1,9 


o-3iiV 


, ,+15.6 

•'•^-o.o 


73 7+19.9 
'3^ '-61.5 


+ 20-7 
"-"-0.0 


113 


K7 


0+1-^ 

"•"-0.0 


+ 0^7 
"•"-0.0 


3+1.7 




24 9 + 12-9 
2^-^-10.3 


io.3l*;9 


5.8l|22 


0+2 
"■"-0.0 


42 , + 13.4 
''2^5-i7,9 


„ + 7.0 
8-8-4.4 






36 9+19-9 
-"'•9-14.1 


2 1 + **^1 
2^^-0.0 


3.6lt;l 


0+2-6 
"-"-0.0 


36.2 + 21-S 


1 7+9.0 
^■'-0.0 


0+58 
"•"-0.0 


0+7^9 
"■"-0.0 


10.9+^138 


a «+10.8 


114 


F9 


4+^-3 
-0.0 


o + o^ = 

"•"-0.0 


2.8i;;2 


o-iiH 


34 7+17-7 
^ -12.5 


4 2+3^2 
^■'^-2.8 


1 9+66 
9-0.6 


0+2^3 
"■"-0.0 


51.S + 15^2 


I7.4IIO2I 


115 


MO. 5 


„ ^ + 6.8 
-2.2 


4 + 1^2 
"•*-0.2 


1 5+1^^ 
^■"-0.0 


o-oi;,.?, 


„_ 7+17.3 
6°-7-48.2 




0+'''3 
"•"-0.0 


ooio.o 


25.6lli0:3 


2-6i;;:^9 


117 


K2 




1 4+3^6 


^■911.1 


0+1-6 
"-"-0.0 


5 1+6-4 


10.7l23;l 


0+0'9 
"•"-0.0 


+ 09 
"■"-0.0 


2 1 + 2-0 
2-^-0.0 


2+1-2 
"-2-0.0 






37 9+l''-9 
■^'-9-15.7 


"■3-4.7 


i9.ot9;2 


6 1 + ^-3 
"■^-4.1 


23 6 + l*-9 


n R+2^1 
"■8-0.4 


0+6^8 
"•"-0.0 


2 3 + 3^9 
2^3-2.4 


2 4+2-8 
2-^-1.3 


1 7+2-2 
l-'-0.9 


119 


K7 


0+11-2 

"■"-0.0 


7+23-7 
"-'-0.0 


1 1 + 17-4 
^■^-0.0 


1 q + 31.6 
^■"-0.0 


11 1 + 56.7 
ii.i_2 7 


2 + 25-1 
"-2-0.0 


6+*1^3 
"•6-0.0 


n 3+33.6 
"•3-0.0 


84 n + 58-6 
84-0_2s 7 


0+96-2 
"-"-0.0 


122 


MO 


30 7+3'^ 
■^"•'-30.4 


+ ^-** 

"■"-0.0 


1 4 + 10.5 

^■*-o.o 


o + ^^i 

"■"-0.0 


52.6110-1 


6 4+15-0 


+ 10^0 
"■"-0.0 


"■"-0.0 


1 7+17.3 

'-■'-i.o 


7 1 + 14.6 
'--'-l.l 






„, , + 38.5 


2 + ''-'^ 
"■2-0.6 


, , + 19.7 
"■"-0.0 


22 .^ + 36.9 
'-14.1 


27.8I393-9 


9+11-2 
"-"-0.0 


+ l*^3 
"■"-0.0 


0.011075 


, , 0+13.6 
118-4. 3 


5 7+18-0 

'-0.0 


123 


MO 


42 o+^^-O 
^'^-"-19. 7 


7+3-3 
"-'-0.0 


5 4+2-7 


o+i^i 

"■"-0.0 


41 2 + 11-7 
"^-2-16.8 


4 1 + 6-4 


0+'12 
"•"-0.0 


2 7+3^7 
2^ '-1.4 


4.0l6;8 


0+0-5 
"-"-0.0 


125 


MO 


4 9+13-0 


+ 3-3 
"-"-0.0 


2 3+2^5 
^■•^-1.6 


0+19 
"■"-0.0 


56.9 23.8 


0+5-0 

"-"-0.0 


0+''^2 
"•"-0.0 


"■"-0.0 


30-815^9 


5 0+11-2 
■'-"-1.2 






31 0+10-3 


11 e+i-fs 

•'^•"-12.8 


0+09 
"■"-2.6 


0+2^ = 
"■"-0.0 


29 1 + 9-2 
'^■'■^-13.7 


7 O+^-l 
'-"-5.2 


2 + 1^8 
"•2-0.0 


+ 3^8 
"■"-0.0 


16.0l6^1 


5 1 + 3-2 
"■1-3.1 


127 


M2 


1 ^ + 2.9 
^•'^-0.0 


3 4+1-2 


+ 20 
"■"-0.0 


9 + 19 
"-■^-0.4 


o-oi^o:"o 


+ 3-1 
"-"-0.0 


1 5 + 5^7 

^•■'-o.o 


1 3 + 4-4 


91 5 + 71-4 
j±.a_21.7 


+ ''' 8 
"■"-0.0 






"•2-0.0 


0+0-** 

"•"-0.0 


g+i^'^ 
"■■'-0.0 


+ 1-3 
"-"-0.0 


l4.7ll;3 


1 7+''-2 
!■ '-0.0 


0+23 
"•"-0.0 


0+0-9 
"-"-0.0 


„2 , + 20.0 
82-5-55.0 


+ 6^7 
"■"-0.0 


137 




40 7+10-3 
*"-'-12.0 


n s + 6. 

"•'^-o.o 


9+2 9 
"■"-0.0 


+ 2-3 
"■"-0.0 


39 + 10^8 


5 6+3-4 
"■6-4.2 


5.8+11 
— 3.6 


0+3-2 
"-"-0.0 


+ 2-3 
"-"-0.0 


7 1 + 4-7 
'-^--2.1 


142 


M4 


7+6-3 
"-'-0.0 


5 + 2^4 

"•°-o.o 


2 6 + 2-3 
2-6-1.1 


1 3 + 2-9 
^-3-0.2 


^8-7_33;8 


4 0+6^6 
*"-2.6 


1^9-0.0 


4 5+''-9 
4-5-2.7 


25 7+9-0 
23-'-10.1 


1 + 6-1 
"-1-0.0 


144 




n Q+5.3 
"-"-0.0 


0+2^1 
"•"-0.0 


6 1 + 1-6 
"-^-1.3 


4 4 + 2-6 


46 2 + 6-0 
*6-^-23.0 


"■2lt3 


0+30 
"•"-0.0 


0+3-0 
"■"-0.0 


32.0+13/ 


+ 2-6 
"■"-0.0 






,„ + 91.2 
^-•^-9.5 


22 1 + 33-3 
2''-^-11.6 


10.9±26^9 


o+'ii-o 

"-"-0.0 


7 3 + 21.8 
'-3-7.8 


10.5120^1 


, + 10.9 


0+22^2 
"•"-0.0 


8.3+19,:9 

— 3.5 


5.511173 


146 


M4 


6.8t2-7 


+ 0-6 
"-"-0.0 


4 4 + 1-2 


2 2 + 1-1 
2-2-1.1 


34 1 + '^-^ 
■ -6.7 


0+11 

"■"-0.0 


4 1 + 3-0 


0+1^3 
"•"-0.0 


32.6 + 8-3 


15.612-9 






17.2t6-3 


+ *-2 
"•"-0.0 


+ 3-8 
"-"-0.0 


+ 1-^ 

"-"-0.0 


27 7+**^l 
2'^'-8.0 


0+3-8 
"-"-0.0 


2 2 + 1-5 
2-2-0.0 


0+2 6 
"•"-0.0 


52 q + 19-8 
"2--J-10.6 


+ 8-7 
"■"-0.0 


147 




"•"-0.0 


+ 2^9 
"•"-0.0 


, 3 + 2.8 
3-3-1,3 


1 3 + 1-6 
^-3-0.9 


34 1+18.9 
34-1-10.9 


q + 9.0 
-'-8-4.6 


n 0+3- 5 
"-"-0.6 


7+''^8 
"•'-0.0 


50 7+18-1 
""- '-41.6 


0.0lJ5„4 


148 


K7 


'•*-o.o 


2.5i2-7 


3 5+3-'* 
3-^-0.8 


+ 10 

"■"-0.0 


70 6 + l**-2 
'"■"-49.9 


3 7 + 8.0 
3- '-2.0 


0+^-3 
"-"-0.0 


2 5+51 
2'3-0.0 


5 1 + 8-9 
"-^-0.0 


4.6l?;ii 






, „ + 36.6 

^•'^-o.o 


17.811639 


30.8124-5 


1 9 + 10-7 
"-0.0 


2 0+13-4 
2-O-0.0 


10.9ll6,,2 


5 O + l''-^ 
o.u_3 7 


2 9+21.5 
9-0.0 


24 7+18-4 
^*-'-9.6 


1 + '1^7 
"■^-0.0 


149 


MO 


,7 + 22.2 
•"•"-13.3 


7+2-9 
"-'-0.0 


1-^1^1 


+ 16 
"■"-0.0 


24 1 + 10-8 
'^"■^-10.7 


3-9li:6 


0+2-3 
"•"-0.0 


1 3+26 
^••3-0.0 


30.1III92 


4+8^8 
"■"-0.0 






41 7+24-7 
*^-'-28.1 


+ ^-^ 

"-"-0.0 


4.lll0,3 


4.4+11 
— 3.6 


98 3 + 14.5 
28.3_2i 1 


3 7+13.2 
3-'-0.5 


1 1 + 12-7 


+ 83 
"•"-0.0 


8.816-5 


7 + 8.4 
7-8-1.9 



04108+2910 

04200+2759 

04216+2603 

04303+2240 
04370+2559 
04385+2550 



BPTau 
CITau 



C'WTau 
CoKuTau3 



CoKuTau4 



DFTau 
DHTau 

DKTau 
DLTau 
DMTau 



K7 
K7 



K3 
Ml 



MO. 5 
M2 

MO 
K7 
Ml 



1 + 0-4 
"-1-0.0 

I3.2I2-3 

4 7+1' 5 
*- ' -2.4 

14.31113= 

«-9t^l;| 

„+27.5 

-9.0 

,+9.1 
-0.0 
,+7.3 
-0.0 
, + 2.3 
-7.9 
2+11.7 
'^-O.O 

12 5+10-1 
12-3-0.6 
B+19.7 
-0.0 
9+5.4 
■'-7.9 

2 3+5.6 

2- 3-0.0 
+ 12.4 

-7.1 

8-2l^'| 

3 3 + 6.8 

3- 3-0.5 
2 + 9.6 

-6.0 

5t?^l 
13.418-4 

6+4-7 
-2.7 
30 7+26.0 
38-7_23.6 

11 0+6-9 
ii.u_2.4 

51 2 + 22-1 
51-2-24.7 

6+8-0 
"-"-0.0 

32- 01^2:1 

"tl:? 
,+5.6 
-1.5 

, + 13.6 
-0.0 
42.6128.1 

3 4+ = 'l 
3^*'-0.0 

33 7+39- 1 

33- I —0.0 



18.63 

2.4^ 



3.1] 
16.2^ 



13.23 



2.63 
29.9^ 



14.8^ 



15.2^ 



24.5^ 



10.6^ 



13.0^ 



0.1^ 



4.1^ 



0-lZ 
0.1 + 



0.03 
1.2 + 



0.8"' 



5+1-* 
"-"-0.0 

,+6.5 
-0.0 
+ 2.2 
-0.0 

1 2 + 'l-5 

1-8-0.7 
1 + 16.8 
"• -0.0 
n+1.7 
-0.0 
2 + 2.3 
-0.0 
8 + 1-3 
-0.5 

0+2-2 
"-"-0.0 

4+1-0 
"• -0.0 

19 2+15^3 
192-0.0 

1 1 + 0.9 

+ 03 
"•"-0.0 

3 + l'5 
"•3-0.0 

1 + 0'7 
lO-O.O 

0+10 

"•"-0.0 
+ 10.2 
-0.0 
+ 2.7 
0.0 

3 + 13 
"•3-0.6 

2 6 + 38 
2^6-0.0 
9 + 8^8 
"•"-6.1 

+ 1.7 
0.0 

3 + 62 
"•3-0.0 

2 + 2^3 
"•^-0.0 

1 8 + 1^8 
1^8-0.8 
, , + 6.0 
^^•"-1.7 

6 + 0^9 
"•6-0.0 

+ 1^1 

"•"-0.0 

0+1^8 
"•"-0.0 

7+8^9 
"•'-0.0 

0+0'6 
"•"-0.0 

2 + 91 
"•2-0.0 



5.3 
0.0^ 



0.4^ 



2 9 + 1^2 
2'"-1.0 
q + lO.S 
-0.0 

-2ll:^ 
,+3.0 
-1.0 

,3 2+31-5 

7+3.6 

-0.0 

7+1-5 
"•'-0.0 

n+19 
-0.0 
„+3.3 
-0.0 

1 4+1--1 

■ -0.7 

+ 2''^2 
"■"-0.0 
+ 1.3 
"-0.0 
,+1.0 

■ -0.0 
, + 10.3 

-0.0 

1 0+1^ = 
lO-O.O 

2 7+1^6 
2^ '-1.4 

3 3+1.9 
3.8_2 7 

n+22 
-0.0 

o+i^-i 

"•"-0.7 

2 5+2^6 
2"-0.0 
18 , + 119 
18^5-8. 6 
,+ 2.5 
-0.0 
4+6.9 
-0.0 
+ 2.1 
8-0.0 
4 + 2-4 
^-1.2 



4.03 
0.4'' 



2.1^ 



2.7^ 



1-03 
1.6+ 



i.o-f 



3.5^ 



1.0^ 



3.33 
4.4+ 



5.43 

3 4 + ^--' 
''-"-3.1 

1 q+1-3 
l--'-1.2 
3 + 2.3 
3-0-1.2 

, + 2.3 
-1.0 
, + 3.6 
-0.0 

2 2+1''' 
22-0. 6 

0+181 
"•"-0.0 



4.73 

1.2 + 



0+0^8 
"•"-0.0 
+9.9 
-0.0 
+0.9 
0.0 
3 8+8.2 
38-7. 2 

0+22 
"■" — 0.0 

1 8+40.3 
'•"-8.2 

+ 0.7 
0.0 
+ 0.4 
0.0 
+ 0^7 
"•"-0.0 
, +11.4 
'•'-0.0 

+ 22.4 
0.0 

+ 0^3 
"•"-0.0 

+ 0^9 
"•"-0.0 
7 5+99 
'•"-0.0 
+ 2.4 
0.0 



0.03 
0.0 + 



0.03 
0.0 + 



18.3 



0.2"^ 



, + 4.6 
-0.0 
n + 10 
-0.0 

13 5 + 91 
133-6.3 
n + 1.4 
-0.0 
+ 5-5 
■ -0.0 
f 0.6 
-0.0 
7+10 2 
"■'-0.0 
3+3.1 
-0.0 

1 5 + 1^6 
l"-0.9 

0+5 6 
"■"-1.3 

3 3 + 13 
33-0.9 

50li'^ 

, 3 + 12.2 
"■3-4.8 

1 4+1^1 

0+''-6 
"•" — 0.0 



0.23 

0.0 + 



0.03 

0.0 + 



0.0^ 



0.3"f 



58 5 + 22^3 
58^5_i7 8 

48 4+"^'' 
,„ „ + 32.2 
58.6I26.4 

i8.9lll;| 

4.4ll6i5 

74 2 + 25-2 
'*-^-17.9 
on 2 + 16.5 
8"-2-33.6 
1+12.3 
-12.6 
,+34.5 
-12.4 



11 2+3-1 
11-2 -4.7 



0.0 



55.13 
47.7+ 



51.8 



30.6l^'5 
9 3 + 20:9 
"•3-0.0 

25 4+5'6 
2"^*-8.4 
. + 10.5 
-12.3 

14 5+12 3 
106-4.3 

23.5IIO35 

63 , + 15-4 
63-3_24.8 

24 1 + 10-5 

44 8 + 10-8 
"*-°-12.1 
, 0+4.0 
"-8-2.6 

50 2 + 19-1 

— 24 5 

16.ll?-t 
, + 6.9 
19.5 

15 1 + 11-0 
l"-l-4,2 

, + 3.8 
7.8 
. + 9.8 
10.3 
,7 3 + 20.6 
l'-3-11.8 
+ 10.1 
12.9 
, + 17.1 
-25.0 
+ 8.5 

13.8 
+ 2.6 
-1.0 
. + 11.4 

3.5 
+ 53.9 
0.0 



37.3 + 



19.2 
29.6 



54.6 
38.0 
18.0 
2.0 
27.6 
46.4 



0.0 



0.. 



+8.0 
-0.0 
7 9+6^2 
'•9-4.1 

0+9^8 
"-" — 0.0 

6 4+17 
"•*-5.1 

2 1 + 181 
21-3. 1 

3 ,+8.5 
"•"-0.0 

3 + ^0 
"•3-0.0 

4 8++^9 
+ 16.4 
-0.0 

3 4+2^6 
3^''-2.6 
, 3+16.8 
'•3-0.0 

+ 1-1 
-0.0 

4 6+6-6 

1 9+191 
1-9-0.0 

+ 3.5 

0.0 
3 7+4 5 
"• '-2.0 
+ 3.0 
-0.0 
3 7+3^5 
3' '-1.4 

0+0^7 
"•"-0.4 
2 + 5.6 
-3.0 

1 + 1^1 
"1-0.2 

7 2+3^8 
'•2-2.2 

+ = ^3 
"•"-0.0 
7+3.9 
-0.0 
+ 5.6 
"-3.1 

5 7+17^2 
"• ' -5.9 

3 1 + 27 
31-1. 9 
,+10.2 

-5.6 
« + 6-6 
-7.4 
o 8+2-8 
0-8-0.0 

3 1 + 2-7 
3-1-2.0 
0+38-7 
"■"-0.0 



1.2"' 



1.8^ 



5.2^^ 



3.73 
11.9"^ 



10.13 
16. 6^1 



3 6+3-1 
3-6-1.6 
0+7-8 
"-0.0 
a+5-3 
-0.0 

0+18-6 
"■"-0.0 
n+2-5 
-0.2 
,+20.5 

-2.9 
„ + 10.7 

-6.3 
2 + 5.7 
-0.0 
7+3.8 
-2.4 
n+9-1 
-0.0 



0.03 
2.8+ 



0.0^ 



1-53 
9.6 + 



1.23 
3.7+ 



0.0"' 



2.613-3 

4 5 + 6-0 
"-"-0.0 

0+0-8 
"-"-0.0 

2 8+3-7 
2-^-1.3 
n+3-4 
-0.0 

"ti:S 

o-4i-^:^o 

14.3 + 



0.03 

0.0+ 



.3 



"-4.3 
9+3-9 
-2.8 

1 + 0-7 
"-1-0.2 

2 1 + 5-1 
2-1-1.6 

0+2-2 
"-"-1.2 
,+3.2 



3.9^ 



1.1^ 



1.1 



o-oio.o 

0+*-9 
"-"-0.0 

0+2-2 
"-"-0.4 

4 3+8-0 
"-3-3.2 

0+1-2 
"-"-0.0 

0+6-0 
"-"-0.0 
2 + 5.4 
2-8-1.4 
2-1 + S^O 



0+0^9 
"■"-0.0 
+ 45-4 
'-0.0 
n+2-5 
-0.0 
7 2+31.8 
'-^ — 20.3 

0+2-8 
"■" — 0.0 
,„ ,+ 17.2 
105-6.0 
n+4-0 
-0.4 
n + 2-1 
-0.0 
n+2-5 
-0.0 
19 9+39-6 
19-9-0.0 
, ,+2.6 
2-1-0.8 
n+13-7 
-0.0 

0+0-5 
"-"-0.0 

1+1-5 
"■ -0.0 
7+7.2 
-5.8 
+0.8 
"-0.0 



0.03 
o.o"' 



2.03 

0.0 + 



o.o"' 



7.0^ 



9.73 
0.0+ 



4.43 

2.6+ 



0.0"! 



9.3^ 



1.4 



0.0 
+ 1.9 
0.9 
+ 16.8 
-0.0 



4 + 7.3 
-2.8 
„+4.5 
-2.3 
q ,+6-4 
9-5-4.3 

1 1++-2 

1- 1-0.6 

2+3-9 
"-2-0.9 

n+1-4 
-0.0 
3+15.2 
-0.0 
, Q + 6.5 

2- 8-0.0 

1 1 + 2-2 
1-1-1.0 

0+22-9 
"-"-7.3 
, 8+1-8 
1-8-0.0 
8 8+7-6 
8-8-4.7 

0+2-1 
"-"-0.0 

8-8104^7 
°-o±?.r 



71 + 16.2 
'"^■■^ — 29.0 

12.911=39 

19-6?-^i.6 
10.3l^'3 

27 1 + 9-7 
2'l-9.5 

44 5+31-4 
**-"-17.3 

0+10 

"•"-0.0 

6.4+;;9 

a 1+5-0 
8-1-1.6 

45 7+13-6 
""-'-20.6 
29 6+37-2 
29-6-0.0 

31 6 + 6-8 
31-6-10.3 
,7 8+13.8 
l'^^8-7.8 

27 4+11^9 
2'*-11.2 



4+9.: 
"-19. 3 

26.0l8-| 

22 5 + 15-4 
22-6-10.1 

6 + 6-9 
6-6-0.6 
53.9126-6 

,+ 11. :5 

-0.3 

i-5l*:^ 

31.316-3 

I8-9IIV 
7, 3 + 18.8 
'1-3-17.4 
, + 17.1 
'-22.4 

22 4 + 9-4 
22-"-8.0 
17 2 + 3-1 
1 '-2-3.7 
24 4+12.3 
2"-*-8.5 
,+ 14.4 
10.8 
,+21.5 
-10.3 
, + 23.9 
'-26.6 
,+11.2 
-0.0 



46.4^ 



1.7"' 



44.6^ 



36.8] 
29.7+ 



43.8^ 
3.2+ 



+ 8-4 
"-"-0.0 

11 4+'^-! 
ii-*-6.4 

16 0+''^l 
16'"-11.1 

4 + 2-6 
-1.4 
n+9-3 
-0.0 



0.4^ 



0.0^ 



3.2 



+ 13.5 



5.5"! 



,+9.8 
-0.0 
8 3+2-8 
8-3-2.3 

11 9+2-7 
ll-9_6.2 

3 9+6-0 
3-9-0.0 
, 3 + 5.3 
13-0. 6 
1 2+24^4 
l'2-O.O 
10.312,? 

18^7lii 
20.3136/ 

i8^5ir.*7 

"•3-0.0 

«t^7:? 

, + 5.3 
-3.2 
9+9.0 
-9.0 

0+1-1 

"-0.0 
r + 5.4 
'•"-3.3 
n + 4-6 
-0.0 
n + 7-3 
-0.0 
+ 9-9 
"-"-0.0 

12.4110/ 

4 2 + 1- + 
"-^-1.1 

7+7.4 
'-4.7 
n + 8-3 
-12.1 

6'8i^:^ 
17.01/4= 



9.03 

3.6"^ 



3.53 
22.9"^ 



0-0"^ 



0.02 
0.0"^ 



10. 7^ 



17.0"^ 



Amorphous olivine and pyroxene combined. 

For each object, the first line corresponds to the 



nponent abund. 



, while the second line corresponds to the cold 1 



nponent abundances, when available. 
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TABLE 7 
Continuation 





(0,1 !,ui) 


^o.l ,,.iu) 


(I 


,1 ,,i„) 


((J.l ,,1,1) 


(1,,'-, ,,,„) 


yl.:, ,,,„) 


(1..5 /,m) 


(1,.5 (ii„) 


((,,() ,,„,) 


(6.0 ,iiii) 


DOTau 


MO 


^■^-0.0 


"•" — 0.0 


1 


"-0.0 


3+*^" 
"■3-0.0 


52.2+2»£2 


4 5 + "^'' 
*"-0.0 


, ,+6.3 
"•3 — 0.0 


1 q + B.2 

'•"-o.o 


19.911"/ 


13 1 + "^" 


DPTau 


MO. 5 


ot- ,+13.7 

-^^■^-15.8 


0+1^^ 
"•"-0.0 


1 


,+1.9 
1-0.0 


"•"-0.0 


^"•^t?l:l 


3+2^* 
"•3_o.o 


13-0.0 


0+13 
"•"-0.0 


2 + 1^^ 
"•2-0.0 


11 0+5^^ 
ll"-0.0 






20 1 + 8.9 


9.8j:^;8 


3 


„ + 3.7 
"-2.1 


7 7+3.8 
'•'-3.8 


"•*l^l:°3 


0+^^* 
"•"-1.6 


-•8t|i 


+ 2^7 
"•"-0.3 


4 + 1^2 
"•"-1.0 


3 1 + 2.3 
■'•1-0.0 


DQTau 


MO 


1 

^■S-1.3 


5 5+i'-'' 
■'■''-1.4 


2 


n + 1.1 
"-0.7 


0+"^" 
"•"-0.0 


20.6l««3i 


6 0+3^1 
"•"-2.0 


"•"lo.o 


1 7+2^8 
'•'-1.0 


02 , + 22.9 


"•"lo.o 








"•1-0.0 





„ + 8.0 
"-0.0 


+ 3^3 
"•"-0.0 


30 s + 36-9 
3"^ = -14.5 


""iao" 


4 6+3^'' 
*"-3.1 


0.01/2^ 


2ii±fo:?, 


29.ll21,;2 


DRTau 


K7 


2,.0+li 


"•"-0.0 


1 


„ + 1.5 
"-0.0 


+ "^^ 
"•"-0.0 


31.6 + |^3 


"•slJ:S 


1 6+2^1 


"•"ti-.S 


20.4 + 9;« 


19.81^^4 


DSTau 


K5 


•'•^■■'-6.3 


"•8-0.7 


2 


+ 2.3 
**-1.7 


2 + 12 
"•^-0.0 


17.4lfi^ 


6 2 + 3^* 
"•2-3.1 


2+l^'' 
"•2-0.0 


+ 11 

"•"-0.0 


•36^4ii;,:» 


0+4-5 
"•"-0.0 






22 + ^*-'' 
''''■"-10.3 


+ 3.5 
"•''-1.8 


2 


7+3.3 
'-1.8 


17 2 + l''^" 
1 '2-9.7 


13 4+10^9 


20 1 + 17.0 
2"^1-10.8 


1+36 
"•1-1.4 


7 0+^ 3 
'•"-4.8 


11 + ''^" 
11. "-4. 9 


6 3 + "" 
"•3-4.1 


F04147+2822 


M4 


^"■''-0.0 


"•"-0.0 


2 


9 + 2.5 
•'-0.0 


o+i^* 

"•"-0.0 


71 2+11^3 
'l''-25.9 


o+i^* 

"•"-0.0 


0+4^" 
"•"-0.0 


1 + 2.8 
"■1-0.0 


3^6ig1 


12-0l|^4 


P04192+2647 




3 4+12-^ 
-0.0 


+ 0.8 
"1-0.0 


3 


g + 1.7 
"-0.8 


+ °-'' 

"•"-0.0 


49 0+1" " 
■^•"-11.6 


ft 7+3.6 
"•'^-4.0 


9+33 
"•''-0.5 


7+*^^ 
"■'-0.1 


21^ = t;,'9' 


12 l+'''2 
121-4.0 






19 7+16.8 
-7.4 


+ 9.0 
"•"-1.6 


14 


, + 15.3 
2-7.1 


3 7+9^1 
3^'-4.5 


16 1+121 
i"i-6.0 


^•3±"3:2 


11 7+16.0 
^•'-7.1 


0+5^3 
"•"-0.0 


19.4112,2 


9 7+8.5 
•'• '-4.8 


F04262+2654 




"•"-0.0 


2 1 + 2-0 
''•1-0.9 





n + 1.2 

"-0.0 


3 9+28 
^•J_2.4 


7+59.3 
"3^'-16.5 


I4.9III73 


+ 3^6 
"•"-0.0 


, 7+9.6 
3^ '-3.1 


19 1 + 121 
191-8. 9 


2 6 + 76 
2"-2.2 


F04297+2246 






0+1^^ 
"•"-0.0 


2 


"-1.5 


3+l^S 
"•3-0.0 


34 + ^^3 


9.5_2 g 


2^3it:*9 


+ 3^2 
"•"-0.0 


44 1+9^8 
""•1-26.1 


9+12'5 
"•9-1.0 






24.7t«l/ 


o , + 15.1 
•''•^-0.0 


0. 


, +32.3 
-0.0 


"•"-0.0 


Q + 66.2 
22-8-0.0 


"•"-0.0 


i37ll^5:l 


5 4+35.6 
"•"-0.0 


9, 3 + 29.8 
213-5. 7 


3 , + 23.3 
3^"-0.0 


F04297+2246A 






^•"-1.8 


3 


,+1.0 
-0.7 


"•"-0.0 


45 9 + 9^" 
^"•"-10. 3 


■'•"-3.1 


3 4 + 3" 
^•"-2.4 


2 9 + 23 
2^'-1.5 


23.9 + ^^1 


2 l + 3^" 
21-0. 






„. , + 15.2 
•^^■^-10.4 


'•*-o.o 


2 


+ 2.1 
1-0.0 


+ *'^2 
"•"-0.0 


20.llllil 


"•2i^o:"o 


"•3-0.0 


o + i^" 

"•"-0.0 


11,21^-3 


29 2+15^'' 
2'2-15.7 


F04570+2520 




2 9 + 11-6 


i+i^' 

"1-0.0 





, +3.9 
1-0.0 


+ 2^'' 
"•"-0.0 


•'"•liil'.i 


^^•3i|:i 


6 6 + "^4 
"•"-3.9 


0+4^7 
"•"-0.0 


20.5111/ 


+ 3^2 
"•"-0.3 






i + *^i 

"■^-8.8 







„ + 9.0 
"-0.0 


■»""±3r?3 


+ 36 
"•"-0.9 


0+^3 
"•"-0.0 


13 7+^^l 
13^ '-13.9 


15 7+^^1 
1"'-12,3 


21 5 + 9^9 
''l''-17.6 


9 + 2^6 
■'•"-7.4 


FMTau 


MO 


^''■"-10.3 




2 


,+2.6 
-1.5 


+ 1^^ 
"•"-0.0 


12 4 + 2^5 
'''•*-9.9 


ft 9+3.2 

32-3. 6 


, , + 6,0 
■'■■'-1.4 


"•"IJ-.o 


l3,ol^;3 


5 + 7-4 






25 + *"5^** 

''"■"-o.o 




3 


, + 4.5 
"-0.0 


+ 17.8 
12-0.0 




i^ii^o 


1 3+17.3 

'•"■-o.o 


5 4+14^0 
"•"-0.0 


28-31// 


25-41// 


FNTau 


M5 


11 3+^-'^ 
'-'-■■^-2.0 




2 


7+2.0 
'-1.0 


1 + 2^0 
"1-0.2 


^7.lil4"4 


10.911;? 


1 + 39 
"1-0.0 


3 2 + 3^2 
"•''-2.3 


20.ll?;^ 


O+l"! " 
"•"-0.0 






^"■"-0.0 


9 , + 29.6 

■'••'-0.0 


12 


g + 37.1 

"-0.0 


0+22.9 
"•"-0.0 


19 + 68-3 
i-'"-0-0 


, ,+ 16.2 
13-0. 


2 9 + 143 
•"-0.0 


1 K + 36.4 
1^3-0.0 


1 ft ft + 56. 1 
133_o.o 


,3 , + 102.6 
23^"-0.0 


FOTau 


M2 




o+*^2 

"•"-0.0 


2 


r + 6.5 
"-0.0 


2 + 33 
"•2-0.0 


20 6+l"^'' 


, , + 10.5 
"•"-2.3 


o + 5^* 

"•"-0.0 


, ,+ 13.5 
"•"-0.0 


,2 ,+23.4 
"2^1-18.6 


o+i^'^i 

"•"-0.0 


FPTau 


M4 


25.8 + 9;| 


, 7+2.6 
3^ -2.4 







0+2^1 
"•"-0.0 


5.5 + 11,3 
— 3.6 


0+3^2 
"•"-0.2 


+ 1^^ 
"•"-0.0 


+ 11 
"•"-0.0 


64 1 + 2S-1 
""•1-20.7 


o.ol// 






8,6+|;t 


'•''-3.7 


2 




10 1+6.6 
l"l-0.0 


(! + 5.5 


,„ ,+ 11.9 
1""-10.1 


9 ft + 1^3 
23-0. 


ft 5 + 9.4 
"•"-0.0 


,, , + 18.2 
3-'-'-9.2 


5 0+4^6 
"•"-0.0 


FQTau 


M2 




1 o+^^i 

i"-0.3 





+ 3.3 
*'-0.2 


r, ., + 7.0 
■'•^-2.7 


49.6 + 32;", 


9 1 + 13.6 
•'•1-4.1 


3 + 3^" 
3"-2.4 


4 8 + '''^3 
•°-3.0 


9.0 + 7^" 
— 3.6 


"•"-0.0 


FSTau 


Ml 


^+10.9 


1 6+3^° 
•"-0.6 





, + 1.4 
-0.0 


1 3 + 1^^ 
^•■'-0.6 


42.6ii,i;*, 


0+2^1 
"•"-0.0 


+ ^^^ 
"•"-0.0 


o + ^^5 

"•"-0.0 


17 + 5^3 
i'"-6.1 


4 6 + 5^2 
"•"-3.2 


FTTau 


c 


on 1 +5.5 
'^"■^-3.9 


1 o + i^* 

l"-0.8 


2 


— 1.5 


0+1" 
"•"-0.0 


i3^9i|:| 


3 9+3^1 
3^-'-1.6 


2 1 + 2.1 
21-1.4 


1 + 16 
"1-0.0 


49 7+15.9 
"•'• '-27.4 


6 4+10.3 
"•"-3.6 








o + i"^" 

"•"-0.0 





„ + 4.4 
"-0.0 


2 9+13^7 
^•^-0.0 


48-8lfg» 


0+23^7 
"•"-0.0 


2 9 + 119 
2^"-0.0 


, + 35.3 
"•"-0.0 


12 + 30 6 
12O-0.0 


1ft 6+15^1 
136-0.0 


FVTau 


K5 




"•"-0.0 


3 


9 + 30 


"•"-0.0 


7+19-4 
22-'-0.0 


12,B + 3^7 


9 + 113 
•'•"-0.0 


+ 3.0 
"•"-0.0 


43. 2112^9 


6 9 + 16-1 
"••'-0.0 


FXTau 


Ml 


28.2 + 11 


^•'-0.0 


3 


4+1.8 


0+"^3 
"•"-0.0 


^•^t;"3:l 


6.3 + ^^? 
— 2.5 


+ 21 
"•"-0.0 


o+i^" 

"•"-0.0 


4 + 3-6'' 
"•"-0-0 


ft + 26 
30-1. 8 








"•" — 0.0 


3 


q+8.7 
''-0.0 


4 6+1^ 3 
*"-0.0 


25.3±1<=,<' 


2 + 2^7 
"•2-0.0 


2+*2 
"•2-0.0 


, 7+10.7 
I'-O.O 


,, ,+ 16-6 
dci.J_14 1 


q 6+219 
96-0.0 


FZTau 


MO 


24 8+''-^ 
''4»_10.9 


^•''-o.o 


1 


q+1-0 
"-0.9 


10 7+2-6 
1"- '-2.6 


19.1±1°4« 




1 ,i+2-2 
l-°-0.5 


0+2^8 
"■"-0.0 


25.ol^,;t 


ll.llt" 


GGTau 


MO 


''■'^-3.7 


"•'^-o.o 


1 


6+1-4 
"-0.0 


o+°^8 

"O-O.O 


^•2lli.l 


0+21 
"O-O.O 


i^9i?:l 


0+11 
"•"-0.0 


i6.5l|;| 


15^7ll'4 


GHTau 


M2 




•'-0.9 


2 


1 + 1-1 
-1-3 


"•'-0.7 


27 , + 15.0 
^'•3-11. 6 


^•3ll:? 




4 1 + 3^1 
"•1-2.9 


45 4+26.6 
""•"-14.6 


"•"i*o:? 






o 7+4.5 
"■'^-4.5 


''•"-2.4 


4 


9 + 5-4 
•'-4.4 


2 7+2-7 
2-'-2.3 


11.5 + 3^2 


^•i±ti 


1 0+1^3 
10-1.0 


12 6 + l" 3 
l''"-S.O 


32.4ll/4;3 


19 7+18.4 
!■'• '-13.8 


GITau 


K6 


„+10.0 
-'°-''-10.5 


4 + 3" 

"•*-o.o 


1 


4+1.3 


0+"^" 
"•"-0.0 


39.3+1^ 


5 7+2^6 
"• '-2.7 


1 6+2^0 
1"-1.0 


1 0+1^3 
l"-0.0 


"•4i;:l 


ft 4 + 2-9 
3-4-2.1 






le.oliOjS 


'• -3.8 


11 


, +6.0 
•1-5.2 


1 9 + ^^l 
^•"-1.9 


20.511"/ 


+ 21 
"•"-0.0 


0+2^7 
"•"-0.0 


1 + 2.8 
"1-0.0 


33 9 + 14^4 
33^-'-12.7 


8 1 + 9-2 
"1-4.9 


GKTau 


K7 


64.5t^i,3^ 


"•"-0.0 


3 


, +1.5 
-1.3 


"•"-0.0 


18.0l3-j 


3^7t?:^ 


"•"t°o:S 


1 7+0.9 
!• '-0.7 


7 3 + 1^3 
'•3-1.8 


1 ft+0^3 
13-0. 8 


GMAur 


K3 


•"■■'-11.2 


+ 3^2 

"•"-0.0 


3 


+ 2.8 
1-0.0 


0+13 
"•"-0.0 


16.0±=-^ 


12^2l|1 


14.9l|;«2 


0+1^2 
"•"-0.0 


, 9 + 10.8 
'•2-0.0 


15 4+''^5 
l-'"-6.7 






, r 7 + 29.0 
^■'■'-0.0 


"•^tfl-.o 


0. 


4 + 29.4 


,7 , + 33.6 
1 '3-0.0 




, O + 20.4 
38-0. 


+ 3^6 
"•"-0.0 


17 4+17.6 
l'"-0.0 


33.6132/ 


1 3 + 14.1 
13-0.0 


GOTau 


MO 


12 + 6.-°' 
^''■"-7.6 


"•"-0.0 


3 


4+2.8 


1 + 2^0 
"1-0.0 


"3^2ir7:'8 


3 3+"^" 
3^3-4.1 


"•"t?:*2 


5 6 + 7^3 
"•"-4.8 


-•4t«5.3 


0.0l2;"2 


HKTau 


Ml 


1 o + 2^" 

^■"-0.3 


4 + 2^" 

"•*-o.o 


1 


+ 1.3 
"-0.0 


1 + 1^1 

"1-0.0 


37 5 + 6-1 
3'^"-21.5 


7 7+2^3 
'• '-4.7 


8.2l^;", 


+ 3^3 
"•"-0.0 


1ft 6 + 3^1 
13^6_6^.0 


24 8lti=5 






19.5lf/ 


_ , + 10.3 
• -0.0 


1 


, + 8.8 
-0.0 


11 1 + 54.3 
111-2. 5 


13.5 + 33^1 


6 7+33.1 
"•'-0.0 


9 q+5.6 
2^'-0.0 


7 4+64.7 
'• -0.0 


21.5 + // 
— 6.3 


10.61// 


HNTau 


K5 


13.9i|;^ 


+ "^* 
"•"-0.0 





, + 0.7 
"-0.0 


o + ^^i 

"•"-0.0 


45.5liy 


+ 13 
"•"-0.0 


0+13 
"•"-0.0 


+ 1^3 
"•"-0.0 


22 9 + 4 3 
22. a_4g 


1""±3'2 


HOTau 


MO. 5 


10.813.7 


7+1^" 
"• '-0.4 


2 


j, + 2.6 
"-0.6 


0+"^'' 
"•"-0.0 


34.6112/ 


6 2 + 36 
"•2-3.9 


^•^-o'.o 


4 + 2.2 
"•"-0.0 


43 O + ll^'' 
"^•"-17.1 


0+l"^3 
"•"-0.0 






,„ , + 19.0 


o+i^2 

"•"-0.0 


1 


„ + 4.0 
"-0.0 


+ 27 
"•"-0.0 


37 6 + 23-2 
3'^"-11.3 


"•*-o.o 


3+2.0 

"•3-0.0 


4 2 + '' l 
"•2-3.4 


ft 3 + 7.0 
"•3-3.6 


9 3 + 3.2 
•'•3-5.8 


HP Tan 


K3 


57 4+10.9 

'"■^-lo.e 


"•"-0.0 


1 


+ 1-0 
"-0-0 


"•"-0.0 


22-ll3:« 


i^iti-.i! 


+ 1^3 
"•"-0.0 


9 + 13 
"••'-0.0 


l4.3l^;1^ 


3 4 + 1^3 


Haro6-37 


K7 


11.4^13^5 


"•"-0.0 


3 


4+2-2 


3 0+2^3 
3"-l.S 


45.oli?;i 


10.4l|^34 


26t^o:'o 


7 + "^5 
'•"-4.0 


17.21// 


"•"±J:i 


IPTau 


MO 


28.0±^^^ 


"•8-0.3 


2 


+ 1.4 
■'-1.7 


o + "^-'* 

"•"-0.0 


33.ol?:| 


4.0_2 4 


39il:°5 


0+1^2 
"•"-0.0 


18.71// 


8.5l3;2 






17 0+**^l 
1 '-"-6.0 


"•''-0.0 


1 


7+5.7 
'-0.1 


ft 9 + ''^6 
"•2-2.9 


19 3+13^2 
l"3_s.9 


4 2+13^2 
"•2-0.0 


1 ft+9^0 
13-0. 


13 7+3^5 
13. 7 


32 3 + 11^1 
" 6.9 


1 7+12-0 
'• '-1.4 


IQTau 


MO. 5 


"■3il.i 


o + i^* 

"•"-0.0 


2 


"-1.1 


+ "^° 
"•"-0.0 


30.7l3^3 


ft 4+3-6 
"•1-2.4 


9 7 + 3.5 
2'-1.4 


1 6 + 24 
l"-0.9 


23.7i5i33 


"•7li:l 






24.3±«;?, 


o+i^2 

"•"-0.0 





7+2.7 

'-0.4 


1 6 + 31 
16-1.0 


43 o+l" " 
"3. ■'-12. 9 


0+*" 
"•"-0.0 


6 + 31 
0"-0.0 


o + "^5 

"•"-0.0 


i47i«3:?i 


14.ll// 


ISTau 


K7 


"■"-0.0 


+ ''^3 
"•"-0.0 


6 


R + 2-9 
"-1.8 


11 6 + 3^" 
ii"-3.3 


47 8+13-5 
"'•"-lO-l 


19.0l^^3 


5 1 + 71 
"1-4.2 


, + 4.4 
"•■'-0.0 


9 4 + 3^4 
•'^ -3.0 


+ 7^7 
"•"-0.0 






•^•"-0.0 


"•"-0.0 


10 


. + 55.5 
-0.0 


6 R+66.0 
•^•3-0.0 


26 5 + 12-8 
2"^"-13.3 


+ 76.9 
"•"-0.0 


14 0+'"'^* 
'"•"-0.0 


4 7+66.6 
-0.0 


26 9 + 67.9 
2"^'-6.9 


7 6 + 47-0 
'•"-0.0 


LkCal5 


K5 


9+11.0 
°°-"-15.1 


+ 2^" 
"•"-0.0 


2 


7+3.8 
'-0.0 


"•olS:? 


21 6+ll^^ 
21"-6.1 


7 7+2^7 
'• '-2.7 


6 3+3^5 
"•3-3.0 


+ 0.6 
"•"-0.0 


4 K + 2^7 
4^3-0.0 


ft + 3.1 
"•3-0.0 






22 0+11^7 
''''■"-8.4 


9 2 + '5 ^1 
•'•''-6.1 





+ 2.0 
"-0.0 


10.2+7;| 


l6.ol^;^ 


6 1 + 5-6 
"•1-3-6 


o + i^5 

"0-0.4 


6 6 + 63 
"•"-4.9 


1 6 9 + 8. 1 
l"''-6.9 


13 7+9-5 
13- '-7. 


RWAur 


K3 


26 1 + ^^" 
'"'■1-9.6 


3+l^" 
"•■'-0.0 


1 


r + 2.1 

■'-0.0 


+ "^^ 
"•"-0.0 


39.0II3I1 


2 3 + 3^0 
''•^-0-0 


+ 13 
"•"-0.0 


+ 1^5 

"•"-0.0 


12.4l-3;3 


1ft 3 + 6-7 
13-3-6.8 


RYTau 


Gl 


6.5l-^;^ 


s.4ii;; 


1 


, +0.8 
1-0.6 


8.0+J^^ 


69 7+13.3 
"•'•'-10.2 


1 9 + 3-3 
'••'-0.0 


+ 16 
"•"-0.0 


1 + 21 
l"-0.0 


3.3li;3 


0+"^2 
"•"-0.0 


SUAur 


Gl 


20.2 + =,2 


"•"-0.0 


1 


, + 1.0 
3-0.9 


+ "-" 

"•"-0.0 


49.7+l^i; 


1 2+1^3 
12-0.0 


2 4 + 3^3 
''•"-0.0 


"•"-0.0 


3 9 + 19 
3.a_l 7 


91 3 + 31 
21-3_3 6 


UYAur 


K7 




"1-0.0 





+ 1.3 
"-0.0 


o+i^-' 

"•"-0.0 


54 3+ll-'l 
'"'•3-13.1 


0+1^3 

"•"-0.0 


i3t?;«o 


+ 1^3 
"•"-0.0 


14 1 + 3^9 
l"^l-6.2 


^"•"-o.o 


V710Tau 


Ml 


27 7+l'»-ii 
'"•'-19.2 


o+i^" 

"•"-0.0 


2 


„ + 2.6 
''-1. 8 


q+1-45 
"••'-0.6 


39,1 + 17.3 


ft 5 + 6-9 
"•"-5-0 


"•"11:^5 


ft + 6^2 
3"-4.6 


3 7+1^9 
3^'-1.8 


3 3+2.5 
33-1. 3 


V773Tau 


K3 




, + 5.9 
38-0. 


1 


, + 2.3 
■'-0.0 


0+1^3 
"•"-0.0 


40,9ll,-i 


4 3 + 9-4 
"•3-0.0 


+ 1^3 
"•"-0.0 


+ 3^" 
"•"-0.0 


lo.oll;^ 


, 9 , + 5.6 
125-6.3 


V836Tau 


K7 


1 2 + ^-* 
-'•''-0.0 


q + 0-8 
"••'-0.6 


4 


5 + 1.9 
-0.8 


0+"^^ 
"•"-0.0 


32.4+;«;» 


9 l + 4^2 

— 3 3 


9 ft + 3^4 
23-1. 9 


3^« + 2-8 


^•1±|5^6 


+ 3.7 
"•"-0.0 






30.7t«q* 


0.0 + g^„« 


1 


7+9.2 
'-0.0 


10.4+glo^ 


39 ft + 52. 5 
•'^•" — 19.6 


3.9123/ 


i^7i;,^o= 


"•^i^oV 


17.01330/ 


1 5 + 37 
15-1.5 


V955Tau 


K5 


"•oli:? 


"•"tS'.'o 


6 


,+2.5 
3-1.6 


io.7+3;i 


,q , + 15.0 
"^"•3-32.6 


4.6+|;3 


o.o+3;"o 


0.0l2;9 




7 2 + "7 


VYTau 


MO 


''•"-0.0 


0-oti^o' 


3. 


+ 16.1 
3-0.0 


o-ol^o 


18-1 + 8-2 


11-0.1 


"•otiy 


"■" — 0.0 


74 0+l''-l 
'"•"-47.6 


3+13 
"•3-0.0 


ZZTauIRS 


M4.5 


^^•oi|:l 


7+3-° 
"■'-0.0 





7+2.2 
'-0.0 




76.ll-:« 


3+*3 
"•3-0.5 


3.3lt:^9 


"■" — 0.0 


"•oii;? 










0+^-° 
"•"-0.2 


8 


,+6.9 
1-5.0 




39.4111:1 


0+37 
"•"-0.4 


9iil:l 


9 7+9^3 
"•7-6.1 


7^ol|;S 


4 6+*-3 
"•"-2.9 



Amorphous olivine and pyroxene combined. 

For each object, the first line corresponds to the warm component abundances, while the second line corresponds to the cold component abundances, when available. 



On the Evolution of Dust Mineralogy 



TABLE 8 
Continuation 



ID 




on P\i-" 


En ■ 


r<.i 




< )ii i"v I- :■: 


En - 




Sil '/ 


on Pvi /, 








(0,1 i,ni) 


(O.i 


) {O.L ,jin) 


(0.1 i,in) 




(i.r, i,n 






((>,() //iii; 





Upper Scorpius 



PBB2002 J160357.9 


M2 


8 2 + 18-3 
"■■'-2.1 


0+1-'' 

"-"-0.0 


^-1-0.2 


o-oli:«o 


42 9+l*-8 


+ 0-5 


0+"-'' 

"-"-0.0 


o+=*^o 

"-"-0.0 


46 8 + 23-3 
46.S_32 4 


0+10-0 

"-"-0.0 


PBB2002 J160823.2 


K9 


— 0.0 


0-n + n n 

— 0,0 


3.8 + 6-2 

— 0,0 


"-"±;,:i 


27 1 + 2-5 
2'-l-24.1 


5.9+3-0 

— 3,6 


o.o+l-l 

— 0,0 


0.6+10„4 
— 0,0 


60.1 + 33.3 
— 22.7 


0.0+5-1 
— 0.1 


PBB2002 J160900.7 


K9 


17 2 + ^-'' 
^ '-^-4.5 


"-=li:i 


i + 2 
''-"-1.6 


0+0-3 
"-"-0.0 


42 3+14.3 


4 1+3-7 


o-ol^:=o 


o-oli:S 


94 n + 2 .0 
2*--'-16.8 


R + 4-0 
3-0-4.0 


PBB2002 J160959.4 


M4 


^'-'-11.4 




4 8 + 3-S 
*'-°-2.4 


1+2-5 


61 9+16-5 


10.6+12^0 


0-01^:^0 


o-ol^o;'o 


3 1+5.9 


1 2 + 11-7 
^-2-0.2 


PBB2002 J161115.3 


Ml 


20 3+*-2 




1 9+2-5 
l-"-1.4 






3-iil:S 




1-51^4 


o-2i!:l 


5 3+1-3 


PBB2002 J161420.2 


MO 


°-°-1.4 


0.0+°;^ 


1 4+0.9 
- -0.9 


0.0+0-3 


67.7+^7-1 


0.6+0-0 






15.2 + 2-1 


7-oll:l 


PZ99 J160357.6 


K5 


0+3-1 


0+2-» 
"-"-0.0 


2 8+5-1 


0-oi?:? 


6 1 + ^-2 
°- -0.0 


3 8 + 2-0 
3-°-0.0 


2+2-2 

"■2-0.0 


+ 2-6 
"-"-0.0 


«7-il^i;? 




PZ99 J161411.0 


KO 


+ 2-0 

"-"-0.0 


+-i-" 

"-" -0.0 


1 1 + 2-1 


1 4+*-* 

^-*-o.o 


8 8 + 4-7 


3-olE:S 


o+i-* 

"-"-0.0 


o+'i-i 

"-"-0.0 


34-b1«5% 


0+3-0 
"-"-0.0 


ScoPMSSl 


MO. 5 


9-8iI^3 


0+2-^^ 
"-"-0.0 


1 9 + 3-3 
^-■'-1.8 


3-"-1.6 


72 4 + 30.3 
'2-*-61.0 


4-6l^y 


0+0-0 

"-"-1.9 


1 9 + 0-2 


r « + 4-2 
6-616.2 


0+1-5 
"-"-0.0 












Eta Chamaeleontis 












J0843 


M3.4 


,7 2+16.9 
^'-''-ll.O 


1 1+8.3 


4 4+3.3 


1 1+2-8 

^-°-o.o 


44 e+19.9 

44-6-21.9 


r q+7.6 


0+10-2 
"-"-0.0 


1 6+11-6 
l-O-O.O 


13 6+10-2 


0+3-0 
"-" — 0.0 


RECX-5 


M3.8 


36 1 + 23-4 


o+2-£> 

"-"-0.0 




3+14.8 
"-3-0.0 


,+ 15.6 
^2-3_3^5 


23 q + 13-9 


+ 3-0 
"-"-0.0 


8.0+6-0 


7 ,+18.9 
'■1-0.0 


0+10-0 

"-"-0.0 


RECX-9 


M4.4 


42.91231-6 


"-"-0.0 


1 9 + 3-3 


+ 8.3 
^■'5-6.7 


31 «+15.8 
3^-6-24.2 


9 + ^0 

"-■"-o.o 


5 4 + 5-3 
"-^-3.2 


4 6+ 10 -'8 


3 3 + 30 
33-1.1 


8 + 3-3 
0-3-0.0 


RECX-11 


K6.5 


36.9111:1 


1 4+3-4 


3 4+8.4 


2 6+*-2 


„„ + 15.7 
26.8_i4 4 


„ n+10-5 


0+2-7 
"■" — 0.0 


3 4+3-7 
3-*-1.9 


17 6+3^1 
^'■^-6.6 


0+1-0 

"-" — 0.0 



Amorphous olivine and pyroxene combined. 

For each object, the first line corresponds to the warm component abundances, while the second line corresponds to the cold component abundances, when available. 



